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Abstract 

Velocity Model of the Shallow Lunar Crust 

Anthony F. Gangi 
Department of Geophysics 
Texas A§M University 

The travel times of the seismic waves obtained for the Apollo-14 and -16 
Active Seismic Experiments and the Apollo- 16 grenade launches are shown to be 
consistent with a powder- layer model of the shallow lunar crust. The velocity 
variation with depth determined from these data is: v(z) * 110 m/sec 
for z less than 10 meters and v(z) = 250 m/sec for z greater than 10 meters. 

The velocity values found for the 10 meter depth are similar to those found by 
Kovach, et al. (1972). The z^^^ depth dependence for the velocity of the topmost 
layer is that predicted on the basis of a powder layer (Gangi, 1972). The 
Amplitude variation of the direct waves as a function of source- to- receiver 
separation, x, is A(x) * A^x exp (-ax) where 1.5 < n < 2.2 and a = 0.047 
neper/m. 

Velocity- spectra analyses of the direct, surface-reflected, bottom-reflected 
and refracted waves give results that are consistent with the velocity model 
inferred from the traveltime data. 
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Introduction . The data from the Apollo-14 and Apollo- 16 Active Seismic 
Experiments (ASE) as well as the Apollo-16 grenade launches and the Apollo-17 
Lunar Seismic Profiling Experiment (LSPE) were used to study the velocity 
structure of the shallow Lunar crust. 

We found that the powder- layer model - which has a theoretical depth 
dependence for its velocity given by (Gangi, 1972): v(z) = 110 m/sec, 
for z in meters - is consistent with the traveltime data for depths of the order 
of 9 to 11 meters. Beneath this depth, there is a discontinuous increase in 
the velocity to approximately 2S0 m/sec for the Apollo-14 site and to 
approximately 300 m/sec for the Apollo-16 site. The latter value is not as 
accurate as the Apollo-14 site value because of the data quality. The velocity 
of 250 m/sec at a depth of about 10 meters is consistent with the results found 
by Kovach, et al., (1972) for the Apollo- 16 site. However, the velocity jump 
found in this study for that depth is from about 161 m/sec to 250 m/sec rather 
than the 114 m/sec to 250 m/sec found by Kovach, et al., (1972). 

Results obtaine d. The results of the investigations perfoimed under this grant 
have been presented in two publications (copies are attached as i^)pendicies) : 
'•Velocity Structure of the Shallow Lunar Crust”, A. F. Gangi and T. E. Yen, 

The Moon and the Planets , v. 20, 1979, pp. 439-468 which is given in Appendix 1 
and Velocity Deteimination of the Very Shallow Lunar Crust . Tzuhua E. Yen, M.S. 
Thesis, Department of Geophysics, Texas A§M University, August, 1979, 107 pp + 
xiii which is in Appendix 2. 

In these publicatic^s, the velocity model for the shallow lunar crust was 
detemined and refined using the traveltimes and amplitude data from the Apollo- 
14 and -16 ASEs as well as the Apollo-16 grenade launches. Various data- 
processing techniques were used to iiq)rove the quality of the original data: 

1) the original data were deglitched by hand, 2 ) they were frequency-bandpass 


2 . 


filtered, 3) traces with the same shot/receiver separation were stacked and 
4) variation on the velocity-spectra method of Taner and Koehler (1969) were 
performed on direct arrivals, surface -reflected arrivals, bottom- reflected 
arrivals and the refracted arrivals. The velocity- spectra technique is a 
beam-steering or array -focusing method which gives a maximum response for the 
output of a receiving array when the proper time delays are inserted in each 
element of the array. 

Velocity-Spectra Analyses . In the velocity-spectra analyses, it was assumed 
that the velocity in the top layer of the lunar crust varied as 2 }"^^ and the 
purpose of the analyses was to determine: 1) the refer«ice velocity (tak«i as 
the velocity at 1 km depth in Yen’s thesis and as the velocity at 1 meter in 
this report — the difference between the two values is a factor of /TO or 
3.162), 2) the depth to the discontinuous velocity change and 3) the value 
of the velocity at or below the velocity discontinuity. From the velocity 
spectrum for the direct waves, the reference velocity was found to be 101 m/sec 
(for the 1-meter reference depth; for the 1-km reference depth, the value 
becomes 320 m/sec, see p. 84f., ^pendix 2). This value was also obtained 
using the surface-reflected waves (ibid) . From the velocity spectra of the 
reflected waves (i.e., those reflected from the velocity discontinuity) the 
reference velocity was found to be 100 m/sec (318 m/sec) and the depth to the 
reflector was found to be 8.4 meters (see p. 86f., Appendix 2). The velocity 
spectra of the refracted waves gave 9.4 + 0.3 meters for the depth to the 
velocity discontinuity and 230 + 16 m/sec for the velocity of the refractor 
(see p. 90f. , Appendix 2) . These values are to be compared with the values 
109 m/sec (345 m/sec) for the reference velocity, 11 meters for the depth to 
the velocity discontinuity, and 254 m/sec for the underlying refractor 
found using the traveltnue data (see p. 53f,, Appendix 2, especially, Fig. 26, 
p.64). 
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Apollo -17 LSPE Data . We found thnt it was not possible to incorporate the 
Apollo-17 LSPE data into the program to obtain traveltimes for the larger 
distances. This was due to the fact that the signal -to-noise ratio was too 
small to allow these data to be used. The same was true for the .Apollo-16 
greiades. 

fpollo-16 Gre nade Launches . The Apollo-16 grenade launches did provide useable 
data even though only two of the geophones (geophones 1 and 2) had useable 
amplitudes for all three launches and even though the launch times of the 
grenades were not known accurately. We found arrivaltimes at geophone 1 for 
the launching of grenades 2,3 and 4 to be 177, 151 and 121 milliseconds, res- 
P®ctively , despite the fact that all the grenade launches were from the same 
•location (see Table 14, p. 67, Appendix 2). Nevertheless, the time differences 
between the arrivals at geophones 1 and 2 ( \diich were approximately the same 
for all three launches, see Table 14, ibid) were used to test the velocity 
model as was one (low-quality) determination of the arrivaltime at geophone 3. 
This latter reading was used to determine the velocity of the high-velocity 
layer ( the refractor) vdiich underlies the powder layer at about 10 m depth. 
These data showed that a powder- layer depth of 9 meters, a reference velocity 
of 110 m/sec (350 m/sec), and a refractor velocity of 250 m/sec were consistent 
with the data (see p.66f. Appendix 2), 

Deconvolution . The data on the Apollo-14 and the Apollo-16 ASEs were de- 
convolved to try to improve the determinations of the arrivaltimes of both the 
first and later arrivals. This procedure is used to both narrow the waveform 
of the seismic events in time and to decrease the rise time of the onset of 
the pulses. Narrowing the pulses would decrease the overlap and interference 
of different arrivals while decreasing the rise times would allow more accurate 
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determination of the arrival times. We found that there was little or no 
benefit to the procedure because the signal-to-noise ratio of the original 
data were so low. Mule it was possible to decrease the rise time and decrease 
the pulse widths, the concommitant increase in the noise prevented any im- 
provement in the determinations of the arrivaltimes of the various events. 

This could be predicted from the spectral analyses made of the signal traces 
which showed that the signal-to-nolse ratio of the traces was low. 

^ litude Analvses . The asplitudes of the direct waves were reanalyzed 
assuming that there was eiqjonential attenuation of the amplitudes with 
distance due to absorption (or scattering) in the powder layer. It had been 

as !. for the 1/6-th power velocity variation and as x'^ for the 

constant-velocity model. The measured variation with the source-to-receiver 
separation, x, was found to be x'^-Sto x‘2-2 (see Appendix 1. Abstract and 

P. 453 f.) when no ejqjonential attenuation term is included. When the as- 
sumed anplitude variation was taken to be: 

A(x) » A^x'^'expC-ax) 

(where x is the source-to-recelver separation, A, is a reference amplitude, 
n is the exponent that measures the spreading of the wave surface as it 
propagates away from the source and a is the attenuation coefficient) the 
measured parameters, which were obtained by a least -squares fit to the data, 

ere found to be. n-1.46, A^“71.7 and a*0.047 nepers/meter (see p. 68f., 

Appendix 2, especially p. 74). 

Hus large decrease in amplitude with distance due to the "absorption" 
term (i.e., a"0.047 nepers/m) is not consistent with the high Q values (low 
attenuation coefficients) found by Latham, et al. (1970). It is not proposed 
that the exponential term - exp(-ax) - is due to absorption in the powder 
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layer, but rather is due to the change in waveform that can be expected for 
waves propagating in an inhomogeneous medium; in particular, the powder 
layer with its l/6th-power variation in the velocity with depth. 

Scattering and Full -wave Analysis . No progress was made on the scattering or 
the full-wave arialysis part of the proposed program. We encountered difficulties 
in reading the data tapes on our con^juter system which held up progress and 
took more time that anticipated. Also, those parts of the program that were 
completed took more effort that anticipated, leaving no time to work on these 
two tasks. We feel these are inportant tasks and should be conpleted, if not 
now, at least some tin® in the future. 

Errata . Two errors exist in publication in Appendix 1. The date of receipt 
of the manuscript is given as 3 January, 1973 and tliis should read 3 January, 1979 
instead. Also, the acknowledgement that the work was performed under NASA 
Grant NSG-74-17 was deleted from the manuscript. 

Sumnary . The results on the analysis of the Apollo-14 and Apollo-16 ASE data, 
the Apollo-16 Grenades and their launches, and the Apollo-17 LSPE data show 
that the velocity structure of the shallow lunar crust istl) there is a powder 
layer about 10 meters thick vdiich has a velocity variation given by v(z)»110(z)^'^^ 
m/sec for z in meters and 2) there is a discontinuous increase in the velocity 
from about 161 m/sec to about 250 m/sec at a depth of 10 meters. These results 
were obtained by using the traveltimes of direct and refracted waves and by 
using the velocity spectra of direct waves, waves reflected from the surface, 
waves reflected fran the discontinuous velocity jump at 10 meters depth and the 
waves refracted along the velocity discontinuity. 
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VELOCITY STRUCTURE OE THE SHALLOW LUNAR CRUST 

ANTHONY V. GANGI and TZUHUA E. YEN 
Department of Geophysics, Texas A <4 Af University, College Station, Texas, USA 

? 

(Received 3 January* \ 91fl) 

Abstract. The date from the ApoUo-14 and Apollo*16 Active Seismic Experiments have been 
reanalyzed and show that a power-law velocity variation with depth* v(z) a (0<r< 

10 m)* is consistent with both the travel times and amplitudes of the first arrivals for souroe-to- 
geophone separations up to 32 m. The data were improved by removing spurious glitches, by filtering 
and stacking. While this improved the signal-to-noise ratios* it was not possible to measure the arrival 
times or amplitudes of the first arrivals beyond 32 m. The data quality precludes a definitive distinc- 
tion between the power-law velocity variation and the layered-velocity model proposed previously. 
However, the physical evidence that the shallow lunar regolith is made up of fine particles adds weight 
to the 1 /6-power velocity model because this is the variation predicted theoretically for self- 
compacting spheres. 

The 1/6-power law predicts the travel time, t{x), varies with separation, x, as t{x) = 
using a first-order theory, the amplitude, A{x), varies as A{x) = m > 1 ; the layer- 

velocity model predicts t{x) = t^{x/x^) and A{x) = Ao(x/Xo)‘*, respectively. The measured exponents 
for the arrival times were between 0.63 and 0.84 while those for the amplitudes were between - 1.5 
and - 2.2. The large variability in the amplitude exponent is due, in part, to the coarseness with which 
the amplitudes are measured (only five bits are used per amplitude measurement) and the variability in 
geophone sensitivity and thumper-shot strengths. 

A least-squares analysis was devised which uses redundancy in the amplitude data to extract the 
geophone sensitivities, shot strengths and amplitude exponent. The method was used on the Apollo-16 
ASE dau and it indicates there may be as niich as 30 to 40% variation in geophone sensitivities (due 
to siting and coupling effects) and 15 to 20% variability in the thumper-shot strengths. However, 
because of the low signal-to-noise ratios in the data, there is not sufficient accuracy or redundancy in 
the data to allow high confidence in these results. 


1. Introduction 

The first lunar seismograms recorded by the Apollo- 11 seismometers (Latham et al, 
1970a, b) surprised many seismologiiis. Their unusually long durations (see Figure 1) 
gave rise to numerous theoretical speculations. Proposed mechanisms ranged from second- 
ary-ejecta effects (Latham er al., 1970a; Chang et al, 1970;Mukl)amcdzhanov, 1970) to 
scattering of the waves by shallow internal fractures and inhomogeneities (Latham etal, 
1970a, b) or by topographic irregularities (Gold and Soter, 1970). It soon became clear 
that the secondary-ejecta mechanisms were not viable ones because the same long dur- 
ation occurred for seismograms from moonquakes with foci in the lunar interior. 

Early data indicated that the compressional-wave velocity was very low near the lunar 
surface (~0.1 kms‘‘; Utham etal, 1970c, Sutton and Duennebier, 1970) and increased 
to approximately 6 km s'* at a depth of 20 km (Latham et al, 1970d). Latham etal. 
(1970a, b) showed that the variation of the amplitude envelope with time and distance 
was consistent with a diffusive-scattering mechanism provided the Q of the medium was 
greater than 3000. 

The Moon and the Planets 20 (1 979) 439-468. 0165- 0807/79/0204 -0439$04 JO 
Covyrtght © 1979 by D. Reidel Publishing Co.. Dordrecht, Holland, and Boston. U.S.A. 







Fig. 2. 


Seismic ray paths for a linearly increasing velocity variadon with depth and topographic 
irregularities. (I rom Gold and Soter* 1970). 




Gold and Soter (1970) interpreted the ApoMo-U data to imply that the shallow lunar 

depth and, through computer simulation using ray acoustics, they were able to approxi- 

Tv tt^ r *r ''7 be explained 

by scattering of the nearly vertically-incldent waves by topographic irregularities (Figure 

2). They also showed that the seismic amplitudes are greatly enhanced in such a medium. 

so that It required less power to transmit seismic waves than previously believed 

Kovach ««/ (1971) proposed a layered model with a stepwise-increasing velocity van- 

ation based on the data of the Active Seismic Experiment (ASE) at the ApoUo-14 landing 

Z' . I'’'- ,r ■. ( V.) ofKMra fo, . lop tape, of 8.5 m tWcknei! 

and a Kp - 299 m s for an underlying layer (the Fra Mauro formation) of 38 to 76 m 
thickness. A similar model was used to interpret the Apollo-16 ASE data and gave a 

Gangi (1972) proposed a sclf-compacting-powder model which gives a velocity varying 
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vi un rivsTKucruMioi tin hiiaklow lunak . uu.sr 

as (he slxlh root uf .he depth; in .his model ,he velocity «, ,he lunar surface goes to zero. 

s n urn gives a long duration to the signal by scattering from t(»pogrnphic hregu- 
armes. very low correla.hm between horl.ontal and vertical dlsplace.rlln.s. a d.an 
s Igna ctiveh.pe that varier with sourccto-receiver separation and a varying spectrun. Lr 
the signal duration, The.sc effects have been note.l by Latham r, u/. they 

also are explained by the diffusive scattering model (Utham rr al., I‘)70c) atid the 
surlacc-irrcgularity scattering model (Ciold and Solcr, 1^)70), 

Kovach and Watkins (I ‘)7.i) extended and refined the layered model by incorporating 
e": 

hi r.s 1 I r U ^ «•'"'<' Moon cannot yet 

c resolved (i.c. whether it is smooth as depicted or a stepwise increase) but one simple 

12 s^s r.r" ^ T “"y Pt'Wder such as the Apollo 1 1 or 

than^e and the steep velocity-depth gradient (" 2kms-'km-')'’. However it is not 
expected that a self-compacting-powder layer of 5 km thick would exist on the Moon- if 

lOOnt ^’"***’ “ P'°‘’‘**’'y 

Dainty or aL (1974) performed a detailed analysis of the diffusive-scattering mechan- 
ism and compared their theoretical results both with lunar data and seismic-modelling 

' h. " '>•« of ,1,0 ,c.,„rt„g l.‘o, I, 

“f H ^ scatterers at the base of the layer is - 5 km and the Q 

of the medturn ts 5000. The co.esponding values for a frequency ofl.OHz arri4km 
scattenng-Iayer thickness. - 2 km between scatterers and a Q of 5000. The thickness of 

should hoTd* for b d *” appropriate for the lunar crust. A similar analysis 

c fo bodywave scattering by topographic irregularities; in this case, the 

scaooer ng-layer thickness would correspond to the surface area over which the nearlv 

Ize) ‘^""^Pond ‘o the spacing between surface scatterers (of wave-length 

amn^T Experiments of Apollo 14 

othei Ifn Experiment (LSPE) data of Apollo 17 and 

s . They asumed a layered model and assumed that the first arrivals (beyond about 
lOm) were seismic refractions. They found their travel time data were consistent with a 

(2TSmr‘1n fo n**’' ‘he top layer of 4m thickness, 

i ® of 32 m (thickness of 28 m), (3) 495 m s‘* to a 

depth of 390m, (4) 960m s‘‘ to a depth of 1385 m and (5) 4700m s'* for a depth din 
to a le«‘ 18Wm. However, this last velocity is determined from a single source (the 

phones). The shallower slrueture is obtained from the traveltime data resulting from 1 
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eight cxplosivc'-packuge detonations and the LM ascent; all these sources are within Kin 
of tlic Kcophone array, Cooper vt at. (I‘)74) show lliese data can be lined well wiUi a 
continiMms, linearly-increasing velocity with deplli, namely, V- 105 + 77«c (in s ') 
lor 2 in meters. They also state that "Various power law velocity nutdels can be made to 
fit the observed data , . ,” when only the explosive-package and I M ascent data are used. 
It is clear there is still some i|uc.stion regarding the velocity variation with depth in the 
sliallow lunar crust (2 I ktn). Since the .shallow lunar crust severely niodifies the 
received signals, eveir those Ihuit large distances, it is important to know this shallow 
velocity variation well, Thcrcrore. it is worthwhile to rcanaly/e tlie data to determine 
which velocity variation wiUt depth is the most probable. The data from the Apollo- 1 4 
and Apollo-lb ASK’s have been reanalyzed and the results are given below. 


2. Apollo*t4 and >16 ASE Data 

The data used in this analysis are from the astronaut-activated thumper device of the 
Apollo-14 and Apollo-16 ASE's. In both experiments, three geophones were sited on 
the surface in a linear array with 45.72 m (150 ft) spacing between geophones (Lauderdale 
and Eichelman, 1974). The thumper device was fired at 4.57m (15 ft) intervals between 
the ends of the arrays (see Figure .1). Firings (shots) 5, 6. 8. 9. lO. 14. 15 and 16 of the 
Apollo- 14 ASE misfiled and no data are available for them. For the Apollo-l6 ASE. two 
shots were oi.iitted between geophones t and 2; namely, those at the 4,57m spacing 
from the two geophones. 

The signals from the geophones ate sampled every 1.887 ms. corresponding to a 
Nyquist frequency of about 265 Hz. Because of data transmission limitations, a trade-off 
between sampling rate and the number of bits per sample had to be made. The result was 
that only five bits were available for each sample. In order to cover the maximum poss- 
ible dynamic range with only 32 possible binary numbers, the seismic signals were log 
compressed for large signal levels. The correspondence of the binary-data values (0-31) 
and the voltage from geophone 1, Apollo 1 6 is shown in Table 1. The other geophone 
voltages have similar correspondence with the binary data. With only 32 levels possible 
for the geophone output voltage, the resulting traces will have a coarse character. Tliis 
makes it difficult to obtain accurate amplitude information if no processing of filtering 
is performed on the data. Fortunately. It Is possible to process the data to obtain reason- 
ably accurate amplitude values. 

- In order to achieve meaningful results from the analysis, it was necessary to improve 
the original ASE data. Figure 4 shows three representative traces of the raw data from 
the Apollo- 16 ASE. These data are from the tenth thumper shot and the source-to- 
receiver separations arc 50.29m (1 65 ft). 4.57m (15ft) and 41.14m (135ft) for geo- 
phones 1. 2 and 3 respectively. The thumper-firing time is 1.2 s after the beginning of 
the traces. Wliiie a high signal-to noise ratio (S/N) exists for the shortest separation, the 
S/N for the other two traces is so low that it is difficult. If not impossible, to pick the 
first arrivals or to measure their amplitudes. In addition, goophonc I shows severe 


2 1 20 19 lA SHOT NO. 

o°* 

*■* 45.72m ®*“ 


-45.72 m 




o 

6«op I 


'* '7 16 IS 14 13 ^2°^ 

***** X y ^ e 9 4 3 2 

A ******* X * 


4.57m 

-45.72 m 


O 

Stop 2 


45.72 m 




Correspondence of binary 


Table I 

data values (B.D.) with th. »» u 
phone 1. ApoUo iJ) 


- 2.299 
- 1.279 
-0.711S 
“0.3958 
- 0.2202 
-0.1225 
-0.06817 
- 0.03793 
- 0.02110 
-0.001174 
-0.000653 


- 0.00363 
- 0.00202 
- 0.00112 
-0.00047 

- 0.00000 

+ 0.00048 

0.00111 

0.00200 

0.00360 

0.00648 
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'N|||()n> I «iANc;| a-„, 

thumper rp-ig-iu 





' ’ MIOt 10 



J.IU1M, J . .«,*„, «J h,. ,,„ "« 




vinnn Y siKUcruKi * ok tiik shallow lunak triist 

THUMPER flP-lG-10 


44S 




1 SEC. 



I ig. 5. “IVjjlitchiHl** versions of the twees in I't^ture 4. 

anJ (.^)a number ot'ditTcrent criteria were u.scd simultaneously to identil'y and correct 
the bad sample values. The result ol the 'deglitching' process is shown in Figure 5 for the 
same traces shown in Figure 4. While this improved the records considerably, it is clear 
the S/N ratios for the geophone-1 and -.t traces are still too low to allow positive identifl- 
cation of the first arrivals. 

To imporve the S/N and smooth out the traces, the data were bandpass filtered with a 
four-p».le, anti-aliascd, Butterworth filter 1 2dB/octavc slopes at both low and high fre- 
quencies) which had dB frequencies at IO..S llz and M>.2$ II/.. The result, for the sante 
three traces, are shown in Figure b. While this improved the S/N significantly and 
improved the character «rf the traces (contparc Figures and 6 ), the S/N for separations 
larger than *).|4m (.^0 ft) was still low because of the decrease in the direct arrival's 
amplitude. 
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rHUMPER flP-16-10 FILTERED 


£ 

O' 




1 SEC. 




»«»iupdMMiiiered versions < 








thctlsmk! «TI^e band ol 

spectr- n (for geophone >. s|,ot 10 Apollo 

samples) is shown in Figure 7 While onlv i u- ‘ seconds (1024 

265 H.) is shown there. ifrclL^^S !;ttk 1" “> 

signal energy is contained between 10 and 90 Hz with Z’ 

rt lViI'tl Flgur«Ts7nd b"''"""" 

To lurther improve the data, the tracer with ti.» * * 

- .. rr rrar 








Separation ^ ApoUo-14 Shot Nos.* Apollo- 16 Shot Nos.* 

No. Of - - - - — - - . 


(ft) 

(m) 

traces 

GP-1 

GP-2 

GP-3 

GP-1 

GP-2 

GP-3 

0 

0.00 

6 

21 

11 

1 

19 

11 

1 

15 

4.57 

5 

20 

12 

2 


10 

2 

30 

9.14 

7 

19 

13 

3 

18 

9.12 

3 

45 

13.71 

6 

18 

- 

4 

17 

8.13 

4 

60 

18.29 

6 

17 

7 

- 

16 

7.14 

5 

75 

22.86 

4 

— 

- 

- 

15 

6.15 

6 

90 

27.43 

6 

- 

17 

7 

14 

5.16 

7 

105 

32.00 

6 

- 

4.18 


13 

4.17 

8 

120 

36.58 

7 

13 

3.19 

- 

12 

3.18 

9 


* Thumper-shot numbers which had the proper separation from the three geophones. 
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STACKED, FILTERED AND AMPLIFIED 
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aMumption is reasonably consistent with the results found by Kovach et al t 

by Kovach era/. 

gi«n .ho,.,„.g«,ph„„ JdXhf Slh^ntr on^rl^ '” .'I:' 
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STftCKED. flLTtREO RNO PMPUMeO 120 Hi -SO mi) 
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1 sec. 


H.sm 

- jtMmu 

9. UH 


IT 


»e.29« 



22.6G1 


m. 











■>v>AVv(<Y\l^ 




Fig. 9. Stacked* filtered and amplified ASF profile (Apollo-14 and ApoUo-16 signals combined; 

Bandpass: 20 to 50 Hz). 


result is shown in Figure 8. The second trace in Figure 8 is at the same separation as the 
middle traces in Figures 4, 5 and 6 (i.e., the geophone-2 trace for the 10th thumper shot 
of the Apollo-16 ASH). For this trace, the S/N improvement should be better than a 
factor of 2; however, this degree of improvement was not achieved. Nevertheless, improve- 
ments in S/N were achieved for this trace, and for.Jhe-other^races at larger, separations, 
by the stacking technique. " ■■ 

The result of summing the deglitched traces is shown in Figure 9. These signals were 
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Fig. 10. Singleiteophone profile, filtered and amplified. (Geophone 1. Apollo-16 ASF; Bandpass- 

30 to 50 Hz). 

filtered, before summing, with a 4-pole, anti-aliased, bandpass, Butterv/orth filter with 
3dB frequencies at 20 and .50 Hz. Arrival times can be determined with some certainty 

for separations up to 32.00m (105 ft); it is difficult, if not impossible, to pick arrivals 
beyond that distance. 
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Table III 

Travel timet (mUhecondt) 


Separution 
X (m) 

Measured travel times ^ 



Calculated travel timet 

1 

2 

3 

4 

5 

A 

B 

C 

4.57 

55 

53 

56 

52 


51.7 

44.0 

40.1 

9.14 

91 

91 

- 

87 

99 

92.1 

87.9 

80.2 

13.71 

123 

123 

124 

- 

128 

129.1 

131.9 

120.3 

18.29 

151 

149 

152 

- 

155 

164.2 

175.8 

160.4 

22.86 

_ 

- 

_ 


177? 

197.7 

219.8 

200.5 

27.43 

206? 

230 

196? 

229 

199? 

230.1 

•245.0 

240.6 

32.00 

255? 

274? 

264? 

274? 

- 

261.7 

•260.3 

280.7 

1 

0.76 

0.80 

0.74 

0.84 

0.63 

5/6 

- 

1 


590 

430 

630 

340 

1200 

350 

104 

114 


S Times with question marks (?) indicate cifficult time determinations. 
* Traveltimes of the first refracted wave (earliest arrival). 

1 . Measured from ApoDo*14 and -16 stacked data (3-66 Hz). 

2. Measured from ApoUo-14 (only) stacked data (3-66 Hz). 

3. Measured from Apollo-16 (only) stacked data (3-66 Hz). 

4. Measured from Apollo-14, geophone-2 profile (3-^6 Hz). 

5. Measured from Apollo-16, geophone-1 profile (3-66 Hz). 

A. Self-compacting-powder model; r B 14.S7x''*(ms). 

B. Apollo-14 layered model (Kovach and Watkins, 1973). 

C. Apollo-16 layered model (Kovach and Watkins, 1973). 


a log/log plot would be a straight line whose slope, m, would be determined by the 
power-law exponent (m = 1 — n). 

The slopes of least-squares-fitted straight lines are given in Table HI along with the vel- 
ocity Vq which conesponds to the velocity extrapolated to Zq = 1 km. As indicated earlier, 
it is not expected that the powdered layer would extend to 1 km; therefore, is not an 
estimate of the velocity at that depth but is merely a constant used to characterize the 
velocity. The depth Zq = 1 km is chosen only for convenience; the reference depth could 
have been chosen to be 1 m, in which case, the Vq’s Table HI would be multiplied by 
(0.001 =0.3162. While the measured slopes are variable, they are all consistently 

lower than m = 1 — n = 1, the value that would be obtained for the constant velocity 
model. The measured values tend to cluster near the value predicted by the self- 
compacting-powder model; namely, m = 1 — 1/6 = 0.833. 

The variation in the reference velocity, t>o, is much greater than that of the slopes; its 
values vary between 340 and 630ms~‘. The slope of 0.63 and reference velocity of 
1 200 m s~* for the ApoUo-1 6, geophone-1 profile (Column 5, Table HI) are not very accu- 
rate because there only three good data points (the traveltimes at 9.14, 13.71 and 
18.29 m) for determining these values. It gave the least consistent values forn and Vq. In 
computing the least-squares lines, the questionable data were given a weight equal to one- 
quarter that of the high-S/N data. 

Travel times were calculated from the ApoUo-14 and Apollo-16 velocity models given 
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by Kovach and WatWns (1973). These arc shown in columns B and C of Table 111. In 
column A, the travel times for a powder-layer model with Vq - 350 m s"* and n - 1/6 are 
tabulated. This latter model was an average model found from all the cases treated when 
the velocity exponent, ri, was constrained to be 1/6. Overall, there are not large differ- 
ences between the measured travel times and the calculated travel limes using any of the 
models. However, the biggest differences between the Kovach and Watkins models and 
the measured values occur at the small separations, precisely where the S/N ratios are 
highest and where the travel times can be picked with the greatest certainty. Their models 
can be made to fit the close-in data simply by introducing a thin, lower velocity layer at 
the surface. But it should be recalled that they already have low velocities for the top 
layers (104 and 114ms"‘ for Apollo-14 and -16, respectively) which are relatively thin 
(8,5 and 1 2.2 m, respectively). 

The travel time data for the combined Apollo-14 and -16 stacked traces (Column 1, 
Table III) are shown in Figure 1 1 along with the least-squares-fltted line. These data are 
from the deglitched traces which have been bandpass filtered with a fourth-order, 
Butterworth filter having 3dB points at 3 and 66 Hz. it can be seen from Figure 11 that 
the straight line is an excellent fit to the data and that it would be difficult to change the 
slope from its given value (0.76) to 1 .0, the latter value corresponding to the constant- 
velocity model. Equally good fits of data points to straight lines were found for the 
Apollo- 14-only and Apollo- 16-only stacked data. 

5. Amplitudes 

^ The travel times of the first arrivals over the 0-32 m range do not demonstrate a clear dis- 
tinction between the powered-layer and the layered-velocity models. The data accuracy is 



SEPARATION (m) 

Pjg. 11. Log-log plot of the travel times versus separations for the stacked and filtered traces. 
(ApoUo-14 and Apollo-16 ASP signals combined; Bandpass: 3 to 66 Hz; measured slope: m ^ 0.76 and 

reference velocity: Uq = 590ms '). 
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such that either niutlel can be accepted. To try to distinguish between the two models, 
tlie amplitudes of the first arrivals were measured and compared with the expected dis- 
tance variation predicted by the two models. 

Since the thumper shots give piimariiy vertical forces aiui the geophones are vertically 
oriented, the amplitude of the direct />-wave arrival In the layered model should vary as 
the inverse square of the separation - 1 c.. 


A{x) = /lo(Jf/jfo)’. 


(3) 


for small separations (see, for example, White, 1965; p. 215). On the other hand, for a 
power-law velocity model, the amplitude variation with separation is given by (sec 
Appendix A) 


A(x) = 


Sip) Ah 
Ivx dx* 


m 


(4) 


where S{p)dp is tire energy radiated in a bundle of rays having ray parameters lying 
between p ~ dp/2 and p -b dp/2, the ray parameter is given by p = sin 0 { 2 )/v( 2 ), ${ 2 ) is 
the angle between the ray and the vertical (z) direction, v{ 2 ) is the velocity variation with 
depth and t is the travel time for the ray (with ray parameter, p) which returns to the sur- 
face at separation x. For the self-compacting-powder model, the amplitude variation is 
estimated to be (see Appendix A) 


i4(x) = i4o(At/Xo)'^'*‘'"^''*, m > 1 . (5) 

where m is a measure of the source radiation pattern in the power-law-velocity medium. 
To Insure integrability of 


E = r Sip) dp, 
Jo 


where E is the energy radiated by the thumper source, we find m > I (see Appendix A). 
This indicates the amplitude decrease of the direct wave with separation is less in the 
powder-layer model than that in the constant-velocity model. This is consistent with the 
conclusion of Gold and Soter (1970) based on their analysis for a linearly increasing vel- 
ocity witli depth. 

The determination of the amplitude variation with separation for the Apollo-14 and 
•16 ASE data is more difficult than determining the travel time data because of: (I) the 
coarseness of the amplitude sampling, (2) the variability of the thurnper-shot strengths, 
(3) the variability of the geophone sensitivities (primarily due to siting and coupling of 
the geophones) and (4) the low S/N ratio for the larger separations. The coarseness of the 
amplitude data is significantly reduced by the interpolating effect of bandpass filtering. 
The variability due to the shot strengths, the geophone sensitivities and the low S/N ratio 
are reduced by the averaging inherent in stacking oi summing traces (provided the signals 
are sufficiently coherent for a given source/receiver separation). 

On the basis of the measured arrival times (at least for separations less than 22.86 m - 
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Tuhle IV 

Amplitude data (arbitnuy unitfi) 


A. Ijuiulpaxs: 6611/ 


v(ml 14^ It \4 2i\)^ 14 .^t 


l6-|t 16 !6-3t 


16 14 , 16-^11 


4.57 

9.14 

13.71 

IH.29 

22.86 

27.43 

32.00 


4.62 

3.17 

U3 

0.52 


3.91 

2.08 

f 


CL15 

0 , 10 ? 


4.68 

0 

*> 

« 

0 . 1 ? 


4.19 

3.36 

1.13 

0.52 

* 

0.17 

0 . 10 ? 


1.43 

0.81 

0.56 

7 

0.13 

9 


5.82 

0.76 

0.32 

0.45 

0.18? 

7 

9 


5.34 

3.17 

0,54 

0.42 

0.23? 

0.18? 

? 


5.69 

0.82 

0.31 

0.26 

9 


3.95 

U6 

0.40 


0.30 

9 


Slope 


1.55 


1.83 


0.10 

0.19? 


0.07 

9 


2.15 


1.78 


1.77 2.01 


B. Butulpass; 20' 40IIz 


2.01 1.97 


2.04 


.v(m) I4~lt |4-2<!)^ 14-3t 16 - 1 1 16-2(3)^ 16-3t 


4.57 

9.14 

13.71 

18.29 

22.86 

27.43 

32.00 


2.38 

1.75 

0.56 

0*24 

* 


2.17 

0.93 

* 

* 


0.13 

9 


2.65 

1.49 

0.22 

* 

* 

0.05 

* 


2.26 

1.75 

0.56 

0.24 


0,09 

9 


0.66 

0.38 

0.29 

9 

0.06 

9 


2.94 

0.39 

0 . 20 ? 

0.22 

9 

9 

9 


2.48 

1.70 


0.21 

0.11 

0.10 


2.34 

0.44 

0.17 

0.13 

9 


0.05 

9 


Slope 1.63 1.59 


2.34 - 1.87 - 2.07 - 1,98 


2.00 -2.08 


14-1 means Geophone 1. ApoUo-14 ASE» etc. 

14-1 means stacked traces, Apollo*14 ASH. 

Nllsled sTo*r Geophones 1 and 2. etc. 

No shot available. 

Low S/N ratio. 


T coherency of the signals exists so that averaging of the amplitudes 

should be possible by sununing of traces. The measured amplitudes are given in Table IV 
Both the amplitudes for individual traces and for stacked traces are given. Measurements 
were made on data that had been bandpass filtered by anti-aliased, fourth-order Butter- 
worth niters with 3 dB frequencies of 3-66 Hz and 2(M0 Hz. It can be seen that there is 
a great deal of scatter in the data. Some of this is due to the thumper-shot variability and 
some due to geophone siting, but the major part is due to low S/N ratio and the coarse- 
ness of the amplitude data. Straight lines were fitted, by least-squares, through the data 
points (on a log-’.og graph) and the slopes of these lines are included in Table IV. A repre- 
sentative plot of the ampUtude data along with its least-squares-fltted line is shown in 
Figure 12. This represents one of the most complete sets of amplitude data available for a 
sii^lc geophone; namely, geophone 3 for the Apollo-16 ASE. The original traces were 
bandpass-filtered (3 dB frequencies at 3 and 66 Hz) prior to measuring the amplitudes. 

Because of the low S/N ratio at the larger separations, it is not certain that a straight 
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Fig. 1 2. Log^log plot of the ampUtudes versus separations (Single geophonc amplitudes; Gcophone 
3. ApoUo-16 ASF; Bandpass: 3 to 66 Hz; measured slope: - 2.01). 


line (on a log-log plot) is the appropriate fitting function. While all the data are fairly 
well fitted by the line in Figure 1 2 (with a slope equal to - 2.01 ). it is clear that the two 
largest amplitude values (at 4.57 and 9.14 m). which have the best S/N ratios, suggest a 
lower slope. 

The slopes found for all the cases with fairly good data lie between — 1.5 and - 2.1. 
However, the possible errors on these slopes are of the order of ± 0.5. The fact that the 
slopes are more negative than 1 and close to - 2, the slope predicted by a simple flat- 
layer model, does not mean the amplitude data verifies that model. From Equation (5), 
the slope predicted by the powder-layer model would be more positive tlian - 1. How- 
ever, this equation and the theory used to predict a slope of - 2 for the flat-layer model 
are based on simply fying assumptions; namely, that all the sources are of equal strength, 
the geophones are equally coupled to the regolith, there is no attenuation by absorption 
in either model, there is no energy loss by conversion of p-wave energy into s-wave energy 
(for the powder-layer model) and there are no scatters in the lunar regolith. The latter 
three effects would increase the ampliti-de loss wiUi distance so that the predicted slopes 
(- 2 for the flat layer and -(13-m)/12 for the powder layer) should be considered 
upper bounds on the measured ones. The variability of the thumper-shot strengths and of 
the geophones would increase the scatter in the data. 

While the amplitude data do not preclude cither model conclusively (as they would 
have if the measured amplitudes decreased more slowly than inversely with separation), 
they do favor the powder-layer model. All the loss mechanisms lead to a greater decrease 
In amplitude than predicted by the simple analyses of the two models. However, the 
amplitude data do not show a more rapid decrease than that predicted for the 
homogeneousdaycr models proposed by Kovach and Watkins (1973) while the data 
clearly do show a more rapid decrease with separation than that predicted by the simple 
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6. Geophone-CoupUng and Shot-Strength VariabUity 


Phh€^f=s 

= G,S/|a:, - .T/r, 

^«^::rrj=-=rrci.'S 

( 1 , - or 3), and to the strength of a particular shot sav V Thi« ^ i* • ^ ^ 

necessary because quite clearlv L.h „» ^ “V Sj. This normalization is 

Letting G/Sj =Ao, Equation (6) becomes 

- ^o{GilGjXSjlSj)\xi~xjl'^, 

10gi4</ - •0«^0 + log(C,/C;)-flog(S;/5rf) + ;;,lOg|;j,-;^^,. 
or, for convenience in writing, 

«(/ - tto + gi+Sj + mXn, 

(") 


'^'’the'^lT '”*^^’ " ‘°8(Gr/C;) and s, = log (S,!Sj) 


i*(«o»m,g.») -£ I Wv(«0 + Wjf<^+g,-fSy-tfy)*. 


(9) 
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as a function of these parameters. Tlie result (see AppenUix B) is the matrix equation 

• = A* p, (10) 

where a is a vector whose components depend only upi>n the measured amplitudes (tfy), 
the weights (n^j) and the measured separations (.Vi>); A is a square matrix whose com- 
ponents depend only upon the weiglits and the measured separations while p is a vector 

whose components are the unknown parameters: gi S|* . . . .aoand w. (The detailed 

form ot this equation is given in Appendix B.) The solution to this matrix equation is 

P = (II) 

If there is sufficient ledundancy in the data* the matrix will he well conditioned and non- 
singuhrand will have a stable inverse. 

T.ie weights are established from the quality of the data. The weights for the Apollo- 
16 data are shown in Table V. From the table it is seen that only 14 of the thumper 
shots gave useful amplitude data (shots I, II* 13 * 14 and were not useable) and. of 
these 14 * only three (shots 6. 7 and 17 ) give measurable tirst-arrival amplitudes on more 
than one geophone. (Shots 1 2 and 1 7 gave amplitudes of O.S^) and 0.06‘.\ respectively* for 
the 3“66 11/ bandpassed traces on geophone 2 : all other amplitudes are given in Table IV). 
Therefore* only six amplitude measurements (two for each shot) are available to deter- 
mine the six parameters ao. w* *xji.TD*Si 7 (when geophone 2 and shot 7 arc used as the 
reference geophone and sliot* respectively). With the relative geophone sensitivities* the 
constant do 8nd the exponent m set by these data, the remaining relative shot strengths 
will be determined by the assumed amplitude variation (Fquation (6) or ;/» aid the 
measured amplitude. 

Having only six correlative amplitude measurements to determine six unknowns (by 
means of the linear equations (8) or ( 10)) means Utere is little redundancy in the ampli- 
tude data. Nevertheless* the solution of these six equations in the rix unknowns do con- 
stitute a least -squares solution. Hits is because weigliting factors are used in tlu^ equations; 
the weiglits can be interpreted to mean that more than six measurements of equal weight 
were made, some of which were identical measurements (lo!.« same shot location)* and llte 
results combined together to give a single result of greater weight. 

If we use six available correlative amplitude values* the matrix equation becomes 




<Wj^> (»ei,i7> 







m 



(W|p 0 


gi 





ft 



0 


n 

w J 


- (W||7> 

^ m 
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Data weifihu and leparationt; Apollo*l6 ASI' 



Weifthu: Wfj 


Separations: jc/i 


Shot 
No. ij) 

Geophonc No. (0 
1 2 3 

Cleophone No. (0 
1 2 

3 

\ 

0 

0 

0 

91.44 

45.72 

0 

2 

0 

0 

1 

86.87 

41.14 

4.57 

3 

0 

0 

1 

82.30 

36.58 

9.14 

4 

0 

0 

1 

77.72 

32.00 

13.71 

5 

0 

0 

1 

73,15 

27.43 

18.29 


0 

1 

i 

t 

i 

68.58 

22.86 

22.86 

7 

0 

1 

i 

4 

64.05 

18.29 

27.43 

8 

0 

1 

0 

59.44 

13.71 

32.00 

9 

0 

1 

0 

54.86 

9.14 

36.58 

10 

0 

1 

0 

50.29 

4.57 

41.14 

n 

0 

0 

0 

45.72 

0 

45.72 

\2 

0 

1 

0 

36.58 

9.14 

54.86 

13 

0 

0 

0 

32.00 

13.71 

59.44 

14 

1 

4T 

0 

0 

27.43 

18.29 

64.05 

15 

0 

0 

0 

22.86 

22.86 

68.58 

16 

1 

0 

0 

18.29 

27.43 

73.15 

17 

1 

1 

4 

0 

13.7! 

32.00 

77.72 

18 

1 

0 

0 

9.14 

36.58 

82.30 

19 

0 

0 

0 

0 

45.72 

91.44 


where i = I, 2 and 3 and / = (>. 7, 1 7. In terms of assumed values of Wij and tlie measured 
values of and atj. this equation becomes 


3.734' 


’3.50 I0..35 1.00 0.75 1.00 1.2s' 


ao 

11.710 


30.84 2.62 2.39 3.13 3.48 


m 

0.211 


1.00 0 0 1.00 


gi 

- 1.164 


0.75 0..50 0 


gj 

- 1.592 


1.00 0 


St 

- 0.914 


- - - _ 1.2S 


J|7 
m m 


If we solve this matrix equation, the relative geophone sensitivities and relative shot 
strengths are found to be 

GtICi = 0.724; C,/«2 * 1.40, 

SJS', - 0.803; S„tS, = 0.848, 
and the exponent is 

m = 3.57. 
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Unfortunately, these values appear to be unreasonable; tills is not surprising consider- 
ing the lack of redundancy and quality in the amplitude data. The 30 to 40% differences 
in the relative geophone sensitivities are not too unreasonable, but are higher than 
expected. Also, the 15 to 20^^ variations in the shot strengths are possible, but again seem 
large. The value of the exponent (m = - 3.57) is different by almost a factor of two com- 
pared to the values obtained using singie-geophone profiles and stacked profiles (compare 
Table IV). The 30 to 40% differences in geophone sensitivity have no effect on the 
amplitude variation with distance as determine 1 by a single-geophone profile. The 20% 
differences in shot strengths (of shots 6 and 17 relative to shot 7) would not cause 
appreciable differences in the slopes (or exponents) obtained from single-geophone pro- 
files (provided, of course, that these differences are representative of the differences in 
the other shots). It is concluded that the least-squares analysis given above does not give 
reliable values for the parameters (/n.C./Cj.Cs/C,. 5^/5,. However, the method 

is a valid one and the reason for the unreliability in the parameter values is the lack of 
redundancy and quality in the data. 

While the method is not useful for this data set. it is presented in detail because there 
may be other instances where it would give valid results. It provides a rationale for the 
dedgn of seismic experiments which test amplitude variation with separation when vari- 
ability in source strengths and geophone sensitivities is anticipated (as is generally the 
case). 

The same amplitude analysis could not be performed on the Apollo- 14 ASE data 
because there were not correlative amplitude values for geophones 2 and 3 (due to mis- 
fires and poor signal-to-noise ratios). 


7. Summary 

The data from the thumper shots of the Apollo- 14 and Apollo- 16 ASE’s have been 
reanalyzed to test whether the velocity variation in the shallow lunar crust (depths < 1 0 m) 
can be represented by a self-compacting-powder-layer as proposed by Gold and Soter 
(1972) and Gang! (1972) or by constont-velocity layers as proposed by Kovach etd 
(1971,1972,1973). 

Both the travel times and the amplitudes for the first arrivals were remeasured and 
compared with the values predicted by the self-compacting-powder-layer model proposed 
by Gangl (o(z) * t>o(2/*o)*^*) #nd the layered-velocity model proposed by Kovach et a! 
To improve the quality of the data, they were ‘deglitched* to remove spurious values and 
bandpass filtered. Four-pole, anti-aliased Butterwortli filters with bandpasses between 3 
and 66 Hz and 20 and 50 Hz (3dB frequencies) were used to improve the signal-to-noisc 
ratio (S/N). In addition, traces from different thumper shots and with the same source- 
to-geophone spacing were summed together to improve the S/N. While these techniques 
improved the S/N. it still was not possible to measure travel times or amplitudes of the 
first arrivals for separations greater than 32 m. 

While there is variability in the results obtained (see Tabic HI), the travel times for the 
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direct arrival over a separation of 32 m can be fit by the 1/6-power velocity model. The 
measured values of the exponent for an assumed power-law velocity varied between 
approximately 1/3 to 1/7; that is, 0.67 < I - w <0.86 (see Table 111) where « is the 
exponent for the depth variation of the velocity. The best (or average) model for both the 
ApoUo-14 and ApoUo-16 sites is v{z) 350(z/zo)‘'‘ m s'* for «o = I km or »(z) 1 lOr*'*, 

0<z < lt)m. This is fairiy close to the velocity variation, v{z) ^ 190z*'*, predicted by 
Cangi (1972) on the basis of Gassmann’s analysis (1953) and the measured mechanical 
properties of the lunar soil. 

The measured travel times of the first arrivals over the 32 m separation are in reason- 
able ^eement with the values predicted by the layered model (see Table 111). However, 
the biggest percentage de.iations occur at the two shortest distances (4.57 and 9.14m) 
where the S/N is high and the travel times can be measured most accurately. At th es e 
reparations, the measured arrival times, which are accurate within at least one sample 
interval (or i.89ms), differs from those predicted by the layered model by 10 to 15ms. 
TTte corresponding differences for the power-law model is generally less than 2 ms. While 
this indicates that the self-compacting-powder-Iayer model is probably the correct one, 
the quality of the data precludes a definitive distinction between the two models. 

No comparison was made of the measured travel times with those predicted by the 
linear velocity variation used by Gold and Soter (1970), namely, v(z) = p# + az, because 
it was an assumed velocity variation which is not based on any physical mechanism. The 
travel time relationship for this velocity variarion, t = (2/e) sinh*‘ {ax/2Vo), should also fit 
the data to the same accuracy as that of the layered-veloclty model. It, too, would have 
the largest percentage deviations at the shortest distances. 

An analysis of the ampUtudes of the first arrivals was performed to test the models. 
The predicted amplitude variation with separation,* - assuming no amplitude loss due to 
attenuation, scattering or conversion of p-wave energy into s-wave energy - for the layer 
model is * * while that for the 1 /6-power velocity model is m> 1. The 

measured exponent varied from - 1.55 to - 2.34 (see Table IV) with the average value 
near - 2.0. While this result, at first glance, seems to favor the constant velocity model, 
the fact that there will be amplitude loss due to scattering, attentuation and wave-type 
conversion makes this result more consistent vrith the power-law model. However, the 
large errois in the amplitude data - which are more severe than the errors in the arrival 
times - pie-dude a definitive conclusion regarding which is the appropriate velocity model. 

An aUempt was made to eliminate the errors in amplitude, due to variations in the 
geophone sensitivities and shot strengths, by using a least-squares method. The method 
requires that the signals, from individual shots, be detected on two or more geophones. 
Unf-)rtunately, only tliree thumper shots were detected on pairs of geophones, and no 
thumper shots gave detectable first arrivals on all three geophones. Consequently, there 
was too little redundancy in the data to give reUable values for the relative geophone 
sensitivities, relative shot strengths or the exponent for the amplitude variation with 
separation. Only for the Apollo-16 ASE was there sufficient data to perform this analy- 
sis at all, and it indicated that there could be 30 to 40% variability in the gcophone 
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Appendix A. Variation in Amplitude with Distance in a 
Vertically Inhomogeneous Medium 
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p-dp/2 


Fig. A-l. (a) Ray paths for a power-law velocity variation: viz) Detail of the wav^ 

fronts dw, in a ray bundle. 


per unit change in the ray parameter, p = the ray parameter - sin e(z)lv{z), 9(z) the 
angle between the ray and the verUcal z axis and is measured counterclockwise from the z 


axis, and t<z) = the velocity variation v«th depth. 

The intensity for any point along the ray can be expressed as 


A* »P) “ 


2itx' iw' 

At a fixed depth, z, and for a particular ray bundle centered about ray parameter p, (see 
Figure A.lb) . dx’ 


dw dx 

iw(x\p) = cosfldx'; or — » cosfl -- 


When the ray reaches the surface (z “ 0), d(0) = ir (cos fl(0) = — 1) and 

x\p,z-0) s Jt(p). 

The intensity at the surface receiver then becomes 
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I(X,P) = 


Sip) dp 


(A.5) 


2vx dx ‘ 

Note that the intensity is positive since both S[p) and x arc positive but dp/dx is negative 
at the surface. 

The general relationship for vertically inhomogeneous media 
p = dr/dx 

holds good, where r(x) is the traveltime; therefore, we have 


ft. 

^ 2nx dx* 


(A.7) 


It can be verified directly that p = dr/dx for a velocity variation of the form 

r(x) = Vo(2/2o)" (0<M<1) 

by using (1); the fact that 

dr/dx = (dr/dp)/(dx/dp) (A.9) 

and (2) the parametric equations for the traveltime, r(p), and the source/receiver 
separation, x(p), (see, for example, Kaufman, 1953) 

tip) = C„p-‘‘ x(p) = (1 -«)C„p->'", (A.10) 

where C„ is a constant and equal to 

2 \/t! 7 .o r(l/ 2 « + 1/2) 


C„ = 


nvl'" r(l/2«) 


(A.11) 


The problem that remains in determining the intensity is to express S(p) in terms of x. 
The rays from a source at the surface propagate, initially, vertically; therefore, for a 
vertical-force source, most of the energy will be directed along the z-axis with little or no 
energy propagating along the surface. The rays received at the surface near the source 
correspond to large values of the ray parameter p, because 

p = l/t<rT). (A.12) 

where Zt is the turning depth of the ray; i.e., its maximum depth of penetraUon. There- 
fore, if we assume an asymptotic expansion for S(p) of the form (for p ► 1) 

S(p) ~ p'^fl + 1/p + 1/P* + . .*)'. m> i. 

and use the fact that (from Equation A.IO) 


we have 


p ~ x~", 

S(p) ~ x"*". 


(A.14) 


This indicates that the source radiation pattern contributes an increase to the intensity 






- lunar crust 
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If we take advantage of the fact that 

1 d^i 

^ yj ^ v-(«42) 


the intensity variation with separation x becomes 

/(X) ~ ;f-tltn(J-m)l, 

^ ’ (A. 16 ) 

and the amplitude variation is 

> 4 (jf) <N, -m )]/2 

For the selfcompacting-powder model, « = 1/6 and we have 
A(x) = ^(*oX^/xo)-(»-mv. 2 . ^ > j 

Appendix B. The Relative Geophone Sensitivities. Shot Strengths and the AmpBtude 
Variation; Least-Squares Analysis 
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space where its gradient is zero; i.e., ‘ *“ “ niininium in the parameter 

J3p£ - V V ^ ^ + AA<, + g, + s,) = 0; 

(B.5) 


or 
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where the term inside the brackets is a so.,»~ 

product of the two vectors. ’ symmetric matrix obtained by the dyadic 
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The resulting set of equations can be written in matrix form as 

■ = A-p, 


(B.7) 

(B.8) 


1 


mid 


467 


Vl;i.OliTY STKlJCTURIi Ol- THE SHALLOW LUNAR CRUST 
where a is the vector 

* = I Pn«n 

"=11 “'</««( *0 + m^y + 8< + «>). (B.9) 

the square-symmetric matrix A is defined in Equation (B.7) and the vector p is defined in 
Equation (B.2). 

The parameter values (oq, m,gi and sf) are just the components of the vector p and 
they are determined by inverting Equation (B.8) to give 

P (B.10) 

This is the solution to the problem of determining the optimum (in a least-squares 
sense) parameters and it can be seen from Equation (B.IO) that the accuracy of the sol- 
ution depends upon the stability of the inverse of the matrix A and the errors in the 
vector a. These, in turn, depend upon the accuracy of the measurements of the separ- 
ations, Xfj, and the amplitudes, rty, and the values of the weights, wy. The weights them- 
selves are established by the accuracies of uy and JTy. In the analysis of the amplitude 
data from the Apollo-16 ASE, it was assumed that the separations were measured with 
high accuracy; consequently, the weights were established only on the basis of the accu- 
racy of the amplitude measurements. 

From the form of matrix A (Equation (B.7)), it can be seen that it will have a stable 
inverse if there is high redundancy in the data; that is, the sums of the weights over the 
geophones (subscript r) and over the shot strengths (subscript /) as well as the double 
sums (over i and/) have large values. This will occur when the ampUtude from each shot is 
measured accurately at each geophone; that is, aU the Wy = 1 . Unfortunately, this is not 
the case in the Apollo- 14 or Apollo- 16 ASE’s. 


References 


BuUen, K. E.: 1963, An Introduction to the Theory of Seismology, Third Edition, Qunbridge Unive^ 
sity Press, New York, N.Y., Chapter 8. 

Chang, G. K.. Gunther, P., and James, D. B.: 1970. ‘A Secondary Ejecta ExpUnation of a Lunar Seis- 
mogram', your. Geophys. Res. 75, 7426. 

Cooper, M. R., Kovach, R. L., and Watkins, J. S.: 1974, 'Lunar Near-Sut£aoe Sttuctute’, Reviews of 
Geophysics end Space Physics 12. 291. 

Dairtty, A. M. , ToksOz, M. N., Anderson, K. R., Pines, P.J., Nakamura, Y., and Latham, G.: 1974, 
‘Setamic Scattering and Shallow Structure of the Moon In Oceanui Ptocellarum', The Moon 9 1 1. 

.langi, A. F.: 1972, The Lunar Seismogram', The Moon 4, 40 . 

Gassmann, F.: 1951, 'Elastic Waves Through a Packing of Spheres'. Geophysics 16. 673. See also Geo- 
physics 18, 269 (errata) (1953). 

Gold, T. and Soter, S.; 1970, ‘ApoUo-12 Seismic Signal; Indication of a Deep Layer of Powder’ 

I AO in'?! ^ * * 


Kaufman, H.: 1953, 'Velocity Functions in Seismie Prospecting’, Geophysics 18. 289. 


468 


ANTHONY I , HANOI AND TZUHUA !■:, VKn 


Kovach, R. L., W«fldni,J,s„ and Lanilm T ’ IQ 71 n a , 

ur's'^X"»rs“ 

Utham. G. V.. Ewing, M„ Dorman, J., Press I' Sutton r m < 

Lfludcfdale^ VV. aod Eichelmaji W* iq 74 u itaj * 

New York, of Sound TVanmtssion, McGraw-Hill Rook Co.. Inc 

Sutton, G. H. and Duenriebier F K • 1 07n *iri .it, ’ ’ 

Wl.uflT'lP*”'"' te floptnto of ,h. u„, s,rt,„ run, Solwy.,. 

WW.0. 1. E. «„,o ^ ^ ^ 


T 


appendix II 


VELOCITY determination OF 


THE VERY SHALLOW 


LUNAR CRUST 


A Thesis 

by 

T2UHUA EDWARD YEN 


Submitted to the Graduate Colle^ce of 
.. T j, , . Texas A&M Universitv In 
P la fulfillment of the requirement for the degree 

MASTER OF SCIENCE 


August 1979 


Major Subject* Geophysics 



ABSTRACT 


Velocity Determination of tho Very Shallow Lunar Crunt 

(AufjuBt 1979) 

Tzuhua Edward Yen, B.S. 

National Central University} Rep. of China 
Chairman of Advisory Committee « Dr. A.P. Ganpji 

The very shallow lunar crust (for depths less than 15 m) 
at the Apollo 14 and 16 landing sites has been determined to 
be a self-compacting powder layer overlying a homogeneous 
layer, The velocity function of the powder layer is given 
by V(z)=v^(z/z^)l/6 _ v^= 3401Z0 m/seo Is a reference 

velocity at a depth of km. The thickness of the powder 
layer is approximately 10+1 m. The constant velocity in the 
homogeneous layer is approximately 250+20 m/sec with an 
unknown thickness. 

The data of the Active Seismic Experiments of Apollo 14 
and 16 missions are improved by deglitching, filtering, and 
stacking. However, the measurements of the traveltimes and 
the amplitudes of the first arrivals from the stacked 
profiles are difficult for separations beyond 32.0 m, while 
the traveltime determinations of the first arrivals from the 

single geophone profiles are difficult for separations 
beyond 45*72 m. 

The travel tirng^^ terminations of Apollo 16 ASE grenade 
launchings at Geoplio^ 3 are also difficult because of the 

niECEUHVu gMiwi - vA nuMo 



low an.plltudos. Alao. there Ic a „axi»um of 57 msec time 
offset amohA the instants of frenade launchinRs. The 
amplitude measurements for the rrenade launchin«s are 
impossible at Oeophones 2 and 3. Because of the hl«h noise 
levels existing in the grenade data of the Apollo 16 ASE, 

the measurements of both the traveltlmes and amplitudes ’ 
have low accuracy. 

The slope, m, on the log/lop traveltime ?raph is 
related to the exponent, n. of the velocity function. 
V(z)-V^(2/z„)". That is. the slope, m, is 5/6 for the 
one-sixth power variation (n.1/6) while m equals 1 for the 
constant velocity function (n-0). The measured slope lies 
between O.65 and 0.82. The best-fitted reference velocity. 
V„. at the Apollo 14 and 16 landinR sites are 345 and 
357 m/sec. respectively, by forcinp the slope to be 5/6. 

A homogeneous layer with a constant velocity of 254 and 
302 m/sec at Apollo 14 and 16 landing sites, respectively, 
is found underlying the surface powder layer. However, the 
constant velocity. 302 m/sec. at Apollo 16 landing site is 
questionable because of little and poor data. 

An amplitude anaiysis (Gangi and Yen. 1979) indicates 
that the amplitude/separation variation is 

S > 1. for the powder layer model and x’^ for the homogeneous 
and layered model. These variations are based on the 
assumptions that. 1) there is no energy loss either by 
conversion or by attenuation. 2) there are no scatterers in 
the lunar regolith. 3) all the thumper shots are of equal 


etren^'th, and 4) all the geophones are equally oriented and 
coupled. The exponent of the amplitude variation from the 
measured data, for separations up to 32.0 m, lies between 
>1.5 and -2.5. These values seem to favor the homogeneous 
and layered model t however, the predicted values of the 
exponents for these two models should be treated as upper 
bounds for the measured data. We also find that there are 
30 to 40^ differences in the relative geophone sensitivities 
and 15 to 20^ differences in the thumper-shot strengths. 
These variations seem high but are not too unreasonable. 

The coefficient of attenuation which dominates the amplitude 
variation at larger separations is 0.04? /m. 

A velocity -spectrum analysis indicates that the 
signal-to-noise ratio for the Apollo 14 and 16 ASEs varies 
from 0.5 to 1.5. Velocity spectra of the direct and 
reflected waves for the Apollo ASE data suggest that the 
reference velocity of the powder layer. V^. is approximately 
320 m/sec. Prom the velocity spectra of the reflected and 
refracted waves, the thickness of the surface powder layer 
is evaluated to be approximately 9 m and the constant 
velocity in the homogeneous medium is 230 m/sec. These 
values of the parameters are consistent with those vlaues 

found from the traveltime investigations. All of them are 
within 10 %, 
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CHAPTER I . INTRODUCTION 


Before the first astronauts landed on the surface of 
the moon, NASA had tried sendln^r seismometers to it. hut 
without sucoess. The first seismometer for the Passive 
Seismic Experiment (PSE) was installed on the moon hy 
astronauts Aldrin and Armstrong of the Apollo 11 mission on 
July 20th. 1969. The seismic data sent back showed some 
unusual characteristics which are not found on Earth's 
seismic records (Figure 1, from Toksos et al. 1976). These 
unusual characteristics (summarized by Oangi. 1972) 

1) the long duration of the signal. 

2) the variable character of the signals (variable 
durations, variable onset and shape of the 
envelope) . 


3) the lack of correlation between the vertical 
and horizontal displacement components. 

the variation of the spectrum of a signal over 
its duration. 

5) the variation of the near-surface. P-wave velocity 
(measurements from Aj m/sec to 106 m/sec). 
Motivated by those characteristios. various authors 


The citations on the following pages follow the style 
of Geophysics. 


TI»C niN 


Figure 1 


The A1 inoonquake. May 2J, I970, as 
received by the ALSEP 12 station. 

Here X is the N-S component (S positive). 
^ positive), and Z is^ 

vertical (up positive). Seismogram 
starts at I305 UT. Distance between 
traces is 22 d.lgital units. 

(Prom Toksoz et al, 1974) 
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(Gold and Soter, I9?0j Ganpi, 1972) suppested that the 
topmost layer of the lunar crust consists of compacted 
powder. 

Gold and Soter (1970) investigated seismic-wave 
propagation for an assumed linear velocity increase with 
depth for a self-compacting powder layer; i.e., V=V^(l+az) 
with Vq= 395 m/sec, a=1.97 /km, and z in km. At full 
compaction, at whatever depth this may occur, they assumed 
the P-wave velocity may be as high as 6 km/sec. In such a 
situation, there would be a propagation channel in which 
the velocity gradient is steep between the top and the depth 
of full compaction. 

Gangi (1972) proposed a sixth-power velocity model for 
the lunar surface layer based on Hertz's theory of the 
deformation of spheres in contact and Gassmann's (1951) 
determination of the P-wave velocity for dry, hexagonal- 
closed-packed, self-compacting, uniform spheres. This 
velocity variation has a rapid increase of velocity with 
depth near the surface. Rays in such a velocity variation 
come back to the surface at (or near to) 0° from the 
vertical because the velocity at the surface is zero. This 
would explain the lack of correlation between the vertical 
and horizontal components. 

Kovach et al (19?1» 1972, 197^) took a more 
conventional approach in interpreting the data from the 
Apollo 14 and 16 Active Seismic Experiment (ASE) shown in 


k 

Figure 2 (from Kovaoh et al . I971). They aGaumod the lunar 
crust was horizontally layered «ith each layer having a 
constant veloolty. The velocity of each layer Is determined 
by asDumlng there are refraction arrivals which travel, for 
most of their paths, along the interfaces between the layers. 

Latham et al (1970) explained the long duration, the 
variable character of the .signals and the lack of correla- 
tion between vertical and horizontal components of motion by 
postulating a high quality factor (0=1^3000) for a highly- 
fractured or brecclated crustal layer. They proposed that 
the long duration was due to multiple scattering of the waves 
from the fractures and showed that this would give rise to 
signals similar to those shown in Figure 1 if the intrinsic 
Q of the medium was as high as 3000. 

A comparison of the traveltime curves of the first 
arrival using the models of Gold and Soter (linear velocity 
variation), Gangi (power-law velocity variation) and Kovach 
et al (horizontally-layered crust) is shown in Figure 3 
(from Cooper et al. 1974). All three models can be made to 
fit the travel time reasonably well. Also included in 
Figure 3 is the z* velocity variation suggested by Carrier 

U971) to Gang! and which is described in detail in Gangi 

(1972). 

If the waive velocity increases in a stepwise fashion 
with depth, the traveltime curve for the first arrival will 
be composed of linear segments. If the velocity function is 
continuous, the travel time curve will be a smooth continuous 
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curve. The accuracy and spareenesa of the data shown in 

Figure 3 do not preclude any of the models. The question 

Of whioh. is the correct or best-fitting velocity function 

will have to be answered by more accurate determination of 

the travel time data and- by determining which of the models 

explain most of the unusual characteristics of the lunar, 
seismogram. 

Kovach et al (1971-1972. 1974) do not address the 
question of the long duration of the lunar seismogram which 
lasts_for over 50 minutes. Both Gold and Soter (1970) and 
Gangi (1972) suggest that "random walking" of the seismic 
waves, caused by undulations of the lunar surface, would 
explain this phenomenon. Utham et al (I970) suggest the 
long duration is due to the multiple scattering of the waves 
from the fractures in the shallow lunar crust. 


CHAPTER II, PREVIOUS GEOPHYSICAL 


INVESTIGATIONS 


For Apollo 12 Passlvo Seismic Experiment (PSE). Latham 

et al (1970) constructed three lunar seismic-velocity models 

(Fifure 9. from Latham et al. 1,70) based on the measurements 

of the physical properties of lunar samples to estimate the 
velocity-depth variation. 


Model 1 assumed the same variation in elastic parameters 
With pressure (or depth) as measured in the laboratory on a 
breccia sample. Model 111 was based on the laboratory 
measured properties of a homogeneous, igneous, lunar rock. 
Model ll_was an attempt to combine the properties of the 
ipneous-rock and breccia samples to produce a model that 


would have the elastic parameters of a highly fractured 
igneous material. The traveltimes for the first and second 
arrivals, which were assumed to be the direct P and S waves 
respectively, fall in between those predicted for the 


homogeneous 
model II. 


. igneous model 111 and the fractured igneous 
Thus, the correct model for the upper 20 km of 


lunar material in the vicinity of the Apollo 12 landing site 

must have a velocity-depth function that falls between models 

II and 111 but is closer to the model of homogeneous Igneous 
rock (model III), 

Prom the thumper firings of Apollo lU ASE, Kovach et al 
( 1971 ) obtained the travel time curves shown in Figure 2. 

They found velocities of 104 m/sec and 299 m/sec for the 


TRAVELTtME, sec 
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uioiMiMJbt., aeg 

Figure 4 . Firp and second arrivals from the 
LM impact plotted on the traveltime 
di. tance curves based on the Apollo 
11 lunar sample data. 

(From Latham et al » I970) 
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direct and refracted waves, respectively (See Figure 2). 

Kovach et .1 (1972. 1974) assumed that the lunar erusf 
is hcacgeneous and hcrisontally layered with constant velocity 
and thickness in each layer. Within each layer, the velocity 
■say increase slightly with depth without violating the 
assumptions of the interpretation, but the first arrivals 
observed on the surface must have travelled either as direct 
waves or as head waves, that is. waves critically refracted 
along the top of each layer. Thus, each layer generates a 
straight-line segment in the traVeltlme curve. They obtained 
the velocity structures shown in Figures 5 and 6 (from 

Kovach et al. 1972 and 1974) for the Apollo 16 and 17 landing 
sites, respectively. 

A deep layer of dust on the moon may provide a very good 
seismic-wave transmission channel, m order to investigate 
the seismic properties of such a medium. Gold and Soter (I970) 
performed a computer simulation using ray theory for a linear 
velocity/depth variation. V(z).V^(l.as). If the surface of 
the moon were perfectly flat, a ray. leaving the origin 
Initially at an angle e^. would be reflected back at the 
snFle 6^ at each subsequent step and. after n cycles, would 
have moved a distance nx„ in a time nt„. However, the lunar 
surface in the vicinity of the seismic experiment is not flat 
but is gently undulating as is characteristic of mare regions. 

If the undulations are random, the ray. at the i-th 
reflection, will be reflected from the surface at an angle 
«1 = 
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where is the angle of the random slope encountered 
(positive if tilted up away from the origin). After n+1 
reflections, the ray would have travelled from the origin a 
net distance 

2 n 

R = — 51 tan 0. 

" a i*0 i 

Since the have random signs and magnitudes (up to 
some cutoff CXj^^) • tan would occasionally change sign; 
that is, the ray would be reversed in direction. Thus, a ray 
encountering random slopes would have travelled a 
considerably smaller net distance from the origin in a given 
time than would be the case for a perfectly flat surface. 

The introduction of random walking would explain the long- 
duration characteristic of the lunar seismogreuns. 

Dainty et al (1974) interpreted the long duration of the 
seismic record as a result of scattering in a surface layer 
overlying a non-scattering elastic medium. Seismic-model 
experiments (Dainty et al, 1974) were used to demonstrate 
this interpretation. Figure 7 (from Dainty et al, 197^-) 
diagrammatically illustrates the seismic-modelling apparatus 
used to produce the seismograms shown in Figure 8 (from 
Dainty et al, 1974). Two experiments (Dainty et al, 1974) 
are illustrated! 1) propagation across a plate with grooves 
cut half-way through, and 2) propagation along the edge of a 
plate with holes drilled within a skin depth (for Rayleigh 
waves) of the edge. The first experiment produces a 
seismogram of the nature of that shown in Figure 1 . They 
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Block diagram of model seismic experiments 
on scattering, (a) Retangular plate (body 
) experiment, (b) Surface wave 
scatteiing experiment, (c) Lunar seismo- 
gram simulation experiment. 

(Prom Dainty et al, 19?4) 




I msec 

Figure 8 . Production of the scattered envelope. 

(a) Retangular plate 36 cm wide, 24 cm 
across as in Figure 7(a) without scatters. 

(b) Same plate with 0.6-cm-diameter holes, 
0.56 holes/cm . The P-wave velocity is 
5*6 mm/jLlsj a Poisson's ratio of O .25 is 
assumed. (From Dainty et al, 1974) 





found that most surface-wave energy scattered into body 
waves and the surface wave disappeared • 

Oangi (1972) proposed a 1/6-power velocity-depth 
variation (i.e., V^az^^^) for the lunar surface layer. 

This variation is based on Hertz's theory of the deformation 
of spheres in contact and Gassmann's (I95I) determination 
of the P-wave velocity for dry, hexagonal-closed-packed, 
self-compacting, uniform spheres. Gassmann found 
Vhp . 

Where is a constant. 

By using the values obtained from lunar samples by 
Kanamori et al (1970, 1971). Gangi (1972) found 

600 (z/Zo)^/^ 

where 2^»1 km. This depth dependence of the vertical 
P-wave velocity (Figure 9i from Gangi, I972) gives rise 
to a very rapidly varying velocity with depth. This 
velocity dependence is consistent with the measured velocity 
dependence on pressure for lunar samples (Figure 10; from 
Gangi, 1972). 

The ray parameter, p, which is given by Snell's law, 
p=sin 6 (z)/V(z), is a constant along any ray in a 
horizontally homogeneous but vertically inhomogeneous 
elastic medium. Since, for the power-law velocity model, 
the velocity at the surface (z*0) is zero, a finite ray- 
parameter value requires all rays to be normally incident 
to the surface. This would explain the lack of correlation 
of the vertical and horizontal displacement components. 
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Pipure 9* Velocity-depth variation for 
self-compactinp dry sphere 
(Prom Ganp-i, 1972) 



Picure 10. P-wave velocity vs. pressure 

(Lunar samples 12052(35), 12065(68) 
and 10065). Refsi Kanamori e.t al (1971) 
and Kanamori et al (I970). 

(Prom Gangi, I972) 




1 ? 


In addition, rays emanating from a surface source would all 

be directed vertically into the lunar interior at the 
surface . 

Carrier (I97I) ougReated to Uangi that the velocity 
■nay increase as the fourth-root of depth near the surface: 
i.e.. V=V^(z/a^)l. Gangi (I972) and Cooper (1974) found 
that the traveltime curves with equals 0.748 km/sec and 
0.914 km/sec. respectively, and km will fit the 
results from Kovach's model very well. But the values of 
Vq were so chosen to make a good fit. Because the velocity 
at the surface is zero, this model would also require all 
rays to be normally incident to the surface. 
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CHAPTER III. FIELD PROCEDURES 


The Active Seismic Experiments (ASE) are part of the 
Apollo Lunar Surface Experiment Packa^^e (ALSEP) of the 
Apollo 14 and 16 missions. The purpose of the experiments 
(Lauderdale et al, 1976a) is to study the internal 
structure and characteristics of the moon to a depth of 
less than one hundred meters. The purpose of the Apollo 1? 
Lunar Seismic Profiling Experiment (LSPE) (Lauderdale et al, 
1976b) i T to determine these characteristics to a depth of 
about one kilometer. 

Locations and Site Description:j » 

On Feburary 5th, 1971. the Apollo 14 Lunar Module (LM) 
landed at 3*65® S latitudi and 17.48° W longitude (Figure 11 
and Table 1) (Lauderdale et al, 1976c). At about 1100 m 
east of the landing site is the Cone Crater which is 
located on a ridge of Fra Mauro Formation (Lauderdale et al, 
1976c). Cone Crater is a sharp-rimmed crater. The Fra 
Mauro Formation is an extensive blanketlike deposit lying 
on a broad band around the basin and is interpreted as 
ejected from the impact. Detailed studies (Lauderdale et al, 
1976c) by the geologists indicate that the Fra Mauro 
Formation is mainly composed of moderately coherent 
breccias. The main characteristic features of the Fra Mauro 
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AP-17 


Figure 11. 


IJIap of the locations of Apollo 
1^. 16, end 17 landing sites. 


Table 1. Locations of the landing sites and 

the relative distances to each other. 


Apollo 


Location 


3.65 S| 17.48 W 
8.97 Sj 15.51 E 
20.17 Ni 30.77 E 






formation are rid/^e<? 

few to o , ‘ “« « 

evoral tone of motere hlrh in the violnlty of the 

Oite. They aiso eu..eate<. that the ri.,es were 

^or»ea hy fiowayo of material alony the yrounO .uri^ the 

excavation of the basin by meteorites. 

on April 21st. ,972. the Apollo lb L« landed at 8.97» S 

a itude and 15.51° E lonyitude (Piyure 11 and Table 1) 

(Uuderdale et al. 1976c, which is at about 1007 Pm 

southeast of Apollo 16 landing site. The U, landiny site 

a a the western edye of the Descartes Hlyhlands which is 

2 i-ar central hlyhlands encompassed by hiyhland plains 

a jacent mountainous areas of hilly and furrowed 

et al . 

et al. 1976c). This area, which is 

underlain by impact-yenerated breccias, may represent a 

recant Of an ancient lunar surface sculptured by impact of 

- erial. It may also be attributed to iyneous and volcanic 
activity from within the moon. 

On December 12th. 107 '' ♦k., a .. 

,„o ^ landed at 

20*17 N latitude and 30 77 ® r i 

,An, V lonsitude which is at about 

m and 996 km to the northeast of Apollo 16 and 16 
andiny sites, respectively (Fiyure U and Table 1 ) 

(Uuderdale et al. 1976 c). The LM landiny site is on the 
oar Of a deep valley, the Taurus-Littrow Valley, which is 
at the eastern rim of and is radial to the Serenitatis 
asin (Uuderdale et ai. , 976 c). The Taurus-Littrow Valley 
nterpreted as a deep yraben formed by the Serenitatis 
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impact. EJeou around many orators on the valley floor 

consists Of basalt. It shows that the «raben is partly 

filled by lava flow underlying a relatively thick layer of 
unconsolidated material. 

Apparatus t 

The lunar-Surface-Experiment Package (ALSEP) of both 
Apollo lb and 16 consisted of three fieophones deployed in a 
lx.near array, a thumper, a mortar package assembly («PA) 
Which contained four Rrenades. end the AISEP central 
station (Uuderdale et al. 1976a). The Apollo 14 and 16 
Active Seismic Experiment (ASE) data were obuined from an 
astronaut-activated thumper, four mortar-launched Rrenades. 
and four prenade launchinps. The Apollo 17 ALSEP consisted 
Of four geophones deployed in a T-array. eight explosive 
packages (EP). and the ALSEP central sution (Lauderdale 
et al. 1976b). The Apollo 17 Lunar Seismic Profiling 
Experiment (LSPE) data were generated by eight EP-s weighing 
from 57 grams to 2722 grams (Cooper and Kovach. I975). 

astronaut-activated thumper was a 
short staff (Figure 12i from Kovach et al. 1977) which was 
used to detonate small explosive charges. The thumper was 
so mounted that it was perpendicular to the base plate at 
the lower end of the .staff. An arm-switch and an 
initiator-selector switch were located at the upper end of 
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Figure 12. Drawing of the ASE. 

(From Kovach et al, I972) 


tne staff. 


Pressure-awitch in the base plate detected the 


. - aetect< 

the instant of shot initiate «« a 

^ A cable connected the 
thumper tt> the central -4* <4: 

Of the Shot. 

There were twenty-one thumper shots pUnned for the 

polio 14 asE end nineteen thumper shots planned for the 

po lo 16 ASE (Lauderdale et al, 1976a) But th 

fi « m But thumper shots 

5^ 8, 9. 10. 14. 15. and U of Apollo 14 ASE misfired. 
The amount of explosive used for the thumper shots are 
not hnown. The first thumper shots for both Apollo 14 and 
8 s were fired at Oeophone 3. The Thumper/Ceophone-3 
8 stance was then increased in increments of 4.59 m toward 

U d A (Table Z, from 

Lauderdale et al, I976a)i that is +h 

• thumper shots 11 for both 

the Apollo 14 and 16 ASE'c were firart 

fired at Geophone 2, and 


■V — 








Table 2. The ASE geophone/thumper distances. 
(After Lauderdale et al, l9?6a) 


Distance, m, to — 


Shot No. 


Geophone 1 Geophone 2 Geophone 3 I AP-14 I AP-16 


91.^4 

86.72 
82.30 

77.72 
73.15 

68.58 

64.00 

59.44 

54.86 

50.29 

45.72 

41.14 

36.58 

32.00 
27.43 

22.86 

18.29 

13.71 

9.14 

4.57 

0.00 


45.72 

41.14 

36.58 

32.00 

27.43 

22.86 

18.29 

13.71 

9.14 

4.57 

0.00 

4.57 

9.14 

13.71 

18.29 

22.86 

27.43 

32.00 

36.58 

41.14 

45.72 


0.00 

4.57 

9.14 

13.71 

18.29 

22.86 

27.43 
32.00 

36.58 

41.14 

45.72 

50.29 

54.86 

59.44 

64.05 

68.58 

73.15 

77.72 

82.30 

86.87 

91.44 


Note: Shots 5. 6, 8. 9, 10, 14, I5, and 16 

of Apollo 14 ASE misfired. 


^ CVi VO C^oo Os 
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thumper shots 21 and 19 for Apollo 14 and 16 ASE's, 
respectively, were fired at Geophone 1. However, two 
thumper shots were omitted between Geophones 1 and 2 for the 
Apollo 16 ASEj namely, those at a distance of 4.57 m from 
these two geophones. 

T he mortar-package assembly (MPAK The MPA consisted 
of a mortar box, a grenade-launch-tube assembly, and 
connecting cable (Figure 12; from Kovach et al, I972). The 
launch-tube assembly, which contained four grenades, was 
mounted in the mortar box. The mortar box was deployed at 
approximately 45® to the surface so that the grenades would 
be sent to their maximum distances, ranging from I50 m for 
Grenade 4 to I500 m for Grenade 1. The four grenades were 
identical except for the amount of launching explosive, 
which weighed from 10 grams to 42 grams, and the high- 
explosive charges, which weighed from 45 grams to 454 
grams (Table 3 ; from Lauderdale et al, 1976 a). The MFA was 
located about 14 m north-north-east of Geophone 1 (Figure 13; 
from Kovach and Watkins. 1973 ). The grenade firing 
direction was parallel to the geophone array and toward 
Geophone 3. The geophone/source distances (Table 3; from 
Lauderdale et al, 1976 a) ranged from 61 . 87 m to 1017.42 m. 

The Apollo 14 ASE grenades were not fired because a 
study of the photographs and the astronaut’s descriptions 
of the position of the MPA suggested that the back-blast 
might effect other experiments (Lauderdale et al , 1976 a). 

The decision war that the grenades would tiot be fired until 
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all other experiments fall. Three of four ^rencdes of 
Apollo 16 ASE were launched in a sequence of 2-4-3, 

Grenade 1 was not launched because the angle sensor went 
off scale. 

Geophones . The geophones of both ASE's and the LSPE 
were miniature seismometers of the moving-coil, magnet type. 

The coil was the inertial mass suspended by springs in the 
magnetic field. Above the na:tural resonant frequency of 
the geophone (at 7.5 Hz), the output was proportional to 
the ground motion. The three geophones of Apollo 14 and 16 
ASE's were deployed in a linear array at approximately 3 m, 

49 m, and 94 m from the ALSEP central stations and were 
connected to it by cables (Figure 13; from Kovach and 
Watkins, 1973)* The four geophones of the Apollo 17 LSPE 
were deployed in a T-array at approximately I50 m west of 
the LM and were also connected to the ALSEP central 
station by cables (Figure 14; from Cooper et al, 1974). 

The ASE's and LSPE characteristics are given in 
Tables 4 and 5 (from Lauderdale et al. 1976a, b), 
respectively. The output voltages from the amplifiers for 
a 5 nm zero-to-peak signal at 10 Hz are given in Table 6. 

Note that there is a 0.02 V (or 12^) difference between 
Geophone 2 (GP 2) of Apollo 14 ASE and Geophone 3 of 
Apollo 16 ASE. These differences in amplitude sensitivities 
will cause moderate errors when we measure the amplitudes. 

It is necessary to correct the output voltages from the 


28 




UJ 

O 

< 

h 

UJ LU 

“ >5 

q w >- 

zj33 

O CL Q. 
0. X UJ 
< UJ O 



o-« 


O' 


</ 


v> 


Figure 14. Map of Apollo 17 landing site showing position of geophone 
array (G-1 to G-4) and locations of explosive charges 
(EP 4 etc.) deployed during traverses. Location of 
lunar module is shown as IM. (From Cooper et al, 1974) 


Thfi Apollo and 16 ASE charactori sties* 
(from Lauderdale et al, 1976a) 



Signal-to-noise ratio (rms signal 
to rms noise in dB for a 10-nm 

peak-to-peaJi signal at 10 Hz) 33.6 33.1 32.9 38.9 37.9 45.0 







7*02 7.12 7.13 





Table 6. Output volta^^es from the ASE amplifiers 
for a 5 nm zero-to-peak signal at 10 Hz 


Apollo 16 


Receiving inttnn* 
extends to 16S cm 



Notei- III Pull ring l - pulls one pin to stert SAFEIARM slide timer 

(21 Pull ring 2 - «»lng up ring; rotete W counterciwimisei 
pull pin to relHse SAFEIARM plate 

<11 Pull ring j . pull* ^ Pins to tree firing pin and start 
tntrmn MItry timtr 

'fo 1',. Arming I'equoncr for an 
LLiTh cxpl O'.; I VC paokogc. 
(Prom Lauderdale et al, 
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the amplifiers for the differences in the Rains and 
generator constantB# 

me Explosive Packagen (Rp) . An EP (Figure I5, from 
Lauderdale et al, 1976b) wan activated by removing three 
pun Pina. Removal of the firat pull pin activated the 
SAPE/ARM alide timer. Removal of the aeeond pull pin 
released the SAPE/ARM elide from its constrained SAFE 
position. Removal of the third pull pin removed a 

constraint on the firing pin and activated the thermal 
battery timer. 

The eight EP's were placed on the lunar surface by the 
Apollo-17 astronauts at various locations ^'Pigure 14 j from 
Cooper et al, 1974), and were identical except for the 
amount of charges (weighing from 57 grams to 2722 grams. 
Table 7 \ from Cooper and Kovach, I975) and the preset 
runout time of the timers. The farther the distance, 
the larger the amount of explosive used. The geophone/ 
source distances (Table 7; from Cooper and Kovach, I975) 
ranged from 101 m (GP 2-EP 8) to 2870 m (GP 4-EP 1). 

Log Compression ! 

The seismic data from the geophones were sampled 
every I.9 msec and 8.5 msec (which correspond to Nyquist 
frequencies of approximately 265 Hz and 59 Hz) for the ASE’s 
and the ISPE systems, respectively (Lauderdale et al, 
I976 a,b). To cover a large dynamic range with only a few 
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Table ?. The LSPE geophone/EP distances. 

from Cooper and Kovach, I975) 


EP No. 

Explosive 
Weight, g 

Distance, m, to 

-- 

GP-1 

GP-2 

GP-3 

GP-4 

1 

2722 

2855 

2758 

2818 

2870 

2 

0 

113 

327 

425 

371 

366 

J 

j, 

57 

242 

341 

288 

287 

4 

57 

269 

172 

215 

220 

5 

1361 

2230 

2330 

2290 

2320 

0 

454 

1195 

1 1240 

1195 

1095 

7 

Q 

227 

800 

865 

810 

672 

0 

113 

179 

101 

122 

112 


bits per sample, the seismic data were log compressed. 

The analog output of the logarithmic compressor was 
converted into five-bit binary data for the ASE’s, and into 
seven -bit binary data for the LSPE by an analog-to-digital 
convertor and transmitted to the earth through the ALSEP 
communication networks. 

The binary data stored on the tapes had to be expanded 
from the output, of the log compressor to obtain the 

actual seismic data. The relationships between the binary 
levels and the outputs, of the log compressor are 

given in Table 8 (from Lauderdale et al, 1976a). For the 
Apollo 14 and 16 ASE systems. Equations la and lb give the 
relationships between the output, and input, 

voltages of the "log compressor" (Lauderdale et al, 1976a). 

Vin = ± ^(V^ut’Vl^Aa)* (la) 

For binary levels 0 to I3, the sign was used, and I7 to 
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Table 8. ASE log compression 

1 rrom Lauderdale et al» 1976a) 

binary levels for represent! nir 

.‘.“•sssras's St-'” 

Binary level log compressor output. V 
out 

? 0.059060 

2 .216540 

o . 374020 

4 .531500 

c .688980 

^ . 846460 


1.003940 

1.161420 

I.3I8900 

1.476380 

1.633860 

1.791340 

1 . 948820 

2.IO6300 

T- 2.263780 

Linear portion 2.421260 
of compressor 

2.578740 

2.736220 

2.893700 

3.051180 

3.208660 

3.366140 
3.523620 
3. 681100 
3.838580 

3. 996060 

4.153540 

4.311020 

4.468500 

4.625980 

4.783460 

4.940940 


Negative 
input signals 


Linear portion 
of compressor 


Positive 
input signals 
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31 • the sifn was usedt and 


''in • <''out- 2 -'* 20 )/V 3 (lb) 

for binary levels 14, I 5 , and 16. 

The values of Vj, V^, and for the Apollo ASE 
systems are riven in Table 9 (from Lauderdale et al, 1976a). 
For example, if the binary level is 27 for Geophone 1 of the 
Apollo 16 ASE, then the true input voltaro is 

Vin = exp ((4. 31 1020-4. 55779 )/0. 267730) 


* 0.3978 volts 

For the Apollo LSPE system, the logarithmic compressor 
has the function (from Lauderdale et al, 1976 b) 


where 


''out = 

V: voltage 



+b 


— constant which determines the slope of the 
transfer function 

b: is specified by the DC offset of the compressor 
output and the system noise level 
The values of M and b are determined by calibration 
of the system to provide at least accuracy of the data 
referenced to the level of the input signal. 




332.0 V for all geophones. 



CHAPTER IV. DATA ANALYSES AND INTERPRETATIONS 
OP APOLLO 14 AND 16 ASE'S 

The ASE data are obtained- from» 1) an astronaut- 
activated thumper (thirteen thumper shots for the Apollo 14 
ASE and nineteen thumper shots for Apollo I6’s). 2) three 
mortar-launched grenades for the Apollo 16 's, and 3) three 
grenade launchings for the Apollo I6*s. The geophone/ 
thumper-source distances are increased, in increments of 
^*57 m, from 0 m up to 91.44 m (Table 2t page 23). The 
amount of charge used for the thumper shots is not known. 
The three grenades of the Apollo 16 ASE are launched, in_a— 
sequence of 2-4-3, to distances of 900 m, I50 m, and 300 m, 
respectively. Grenade 1, planned for a launch distance of 
1500 m, was not launched because the angle sensor went 
off scale. The geophone/s ource distances ranged from 
61.87 m to 1017.42 m, and the explosive charges weighed 
from 45 grams to 454 grams for the grenades and from 
10 grams to 42 grams for the thrusts of the grenade 
launchings (Table 3 j page 25). The three grenade 
launchings had the same distances to Geophones 1 (14 m), 

2 (50 m), and 3 (95 m). 

Deglitching* 


The lunar seismograms from the Apollo ASE thumper shots 
and grenades show severe "glitches". 


Most are almost 
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uniformly spaced and of uniform amplitude at binary level 27 
but others sre of varying amplitudes and locations in 
Oeophone 1. while extraneous random glitches exist In the 
data for Geophones 2 and 3 for both Apollo 14 and 16 ASE-^. 

Figure 16 shows three representative traces of the raw 
data from the Apollo 14 and 16 ASE thumper shots and 
grenades. These data are from the tenth thumper shot and the 
geophone/source distances are JO. 29 m. 4.57 m. and 41.14 m 
to Oeophones 1 (OP 1), 2. and 3. respectively. The thumper 
firing time is 1.21 seconds after the beginnings of the 
traces. The data show that there are glitches throughout 
the records 1 they are recognised by the fact that they are 
of short duration - generally, only one sample value - 
and have values which are inconsistent with the preceding 
and/or following sample values. Figure 17 shows another 
three traces of raw data. These data are from the second 
g«nade launching of the Apollo 16 ASE. The geophone/source 
distances are 14 m, JO m, and 9J m to Geophones 1, 2, and 3, 
respectively. The grenade launching time is approximately 
0.60 seconds after the beginnings of the traces. Figure I7 
shows there are relatively fewer glitches for these data 
than there are for the ASE thumper shots (see Figure 16) 
and grenades. Furthermore, no glitches are found in the 
grenade launching data for the second and third geophoncs. 

The first data-improving operation performed is to go 
through the data by hand and remove the extraneous values 


Figure 16. Representative raw data which is from 
the tenth thumper shot of Apollo 16 
ASE. The geophone/source distances 
are 50 . 29 m, 4.57 m, and 41.14 m to 
Geophones 1 (GP 1), 2, and 3 , 
respectively. The firing time is 
1.21 seconds after the beginnings 
of the traces. 






AP-16-fiREHADE LAUNCH-2 



ilk — 

o 

r 

ro 


Figure 17. Representative raw data which is from 

launching of Apollo 
10 ASE. The geophone/source distances 

1?’ 1 ° 95 m to Geophones 

mv . 2, and 3, respectively. 

The firing time is 0.64 seconds after 
the beginnings of the traces. 
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and replace then, by valuee Interpolated from neiyhboriny 
values. A deyUtchln* program Is not used in this process 
because. 1) there are relatively few ylltches (oxcludin* 
the regular, periodic ones in OP 1. there are fewer than 
1^). 2) the coarseness of the amplitude valuee precludes 
automatic, computer interpolation, and 3) a criteria is 
used to identify and correct the bad sample values, that is 
we change the value of the binary level of the bad sample 
by one bit and this makes its amplitude consistent with 
the preceding and/or following sample values. The results 
of the “deglitching" process is shown in Figure 18 for 
the same traces shown in Figure 1$. The seismic traces 
are smoother and much clearer than before. 


^plitude Spectrum Analys e s and Pi Iteri np- . 


For the ASE thumper shots, the inconsistency of the 
amplitude spectra is observed. The spectra shown in this 
section are taken from approximately the first one 
second, 512 sample points, of the data after the firing 
time. Figure 1? shows the amplitude spectrum from the 
second geo.ihone of the tenth thumper shot of the Apollo 16 
ASE. The geophone/source distance is 4.57 m. Most of the 
energy of the signal is concentrated between approximately 
10 Hz and 40 Hz. There are two peaks in the amplitude 
spectrum (with amplitudes of approximately 0.6 and 0.5) at 
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THUMPER flP-16-10 



1 SEC 





Figure 18. Representative of deglitched data 
of the ASE. It is from the tenth 
thumper shot of Apollo 16 ASE. 



-shot 10 of Apollo 16 
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about 18 Hz and 33 Hz, respectively. The noise level 
(with amplitude of approximately 0.1) is reached at 
frequencies greater than 90 Hz. Figure 20 shows the 
amplitude spectrum from the first geophone of the 
seventeenth thumper shot of the Apollo 14 ASE. 

The geophone/source distance is 18.29 m. The energy of 
the signal in still concentrated but harmonically 
distributed between approximately 10 and 40 Hz. The 
maximum amplitude (at 18 Hz) is approximately that of 
the maximum amplitude in Figure 19 . 

The amplitude spectra of the ASE grenades and grenade 
launchings show more consistency than those of the ASE 
thumper shots. The amplitude spectrum has the broadest 
frequency band (approximately from 4 to 30 Hz) for the 
shortest geophone/source distance, 61.8? m; namely, the 
signal at Geophone 3 of Grenade 4. In general, the 
frequency band of the signal narrows when the geophone/ 
source distance increases, both the maximum amplitude and 
the frequency band shift toward lower frequencies. The 
longer the distance the seismic waves travel, the more 
high frequency components the moon filters out. For 
example, the amplitude spectrum of the signal from Geophone 
3 for Grenade 3, whose geophone/source distance is 353.26 m, 
has a frequency band from about 7 Hz to 25 Hz and a maximum 
amplitude at about 1? Hz. The frequency band of the signal 
from Geophone 1 for Grenade 2, which has the largest 
geophone/source distance (1017.42 m), is from about 5 Hz to 
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20 Hi. and Ito maximum amplitude Is at about 11 Ht 

To improve the S/N ratio, the ASE thumper data are 
^ndpaaa mtered by a Pour-pole, antl-allaaed, Butterworth 
ter with -3 dB frequencies at about 3 Hz and 66 Hz 
( .01 f^ and 0.25 f^, respectively, where f„^265 Hz is the 
yqulet frequency). The ASE grenade and grenade launching 
data are bandpass filtered by a similar Butterworth filter 
with -3 dB frequencies at 3 Hz and 4o Hz (0.01 f and 
0.15 f„. respectively). The Butterworth filter (Oppenheim 
and Schafer, I975) has the properties. 1) the amplitude 
response is maximally flat in the passband and 2) the 
approximation is monotonio in the passband and stopband, 
the amplitude, A(n), and the phase, 9(n), functions of a 
2 -point version of a Butterworth bandpass filter are. 

A(n) s ((1+B^)(l+C^))“i 

6(n) - tan (/2C/(l-c^) )-tan"^ (/2B/(1-B^)) 

B=tan (njT/2^'‘h/tan in 
C=tan (7TV2fL)/tan inlT/Z^^h 

"the high and low outoff 
frequencies, respectively. 

0 , 1 , 2 2^-1 

These improve the S/N ratio significantly. 

Representative results are shown in Figure 21 (compare with 

Figure 18). The traces are smoother, and the high frequency 

noise and the low frequency characteristic have been 
eliminated. 
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where 
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WliUe a hiph S/N ratio exiats for tho shortest ASE 
thumper distanoe (4.57 m. OP 2 , Figure 21 ). the S/N ratios 
for the other two traces in Fi,;ure 21 are so low that it 
IS difficult, if not impossible, to pick the first arrivals 
or to measure their amplitudes. The low S/N ratio at 
larger distances is caused by the decrease of the first- 
arrival amplitudes. Kovach et al ( 1972 ) also indicate 
that the onsets become emergent when the separations are 
larger than 15 m and the determination of the signal onsets 
from the background noise is difficult for the Apollo 14 
and 16 ASE thumper shots for separations greater than 45 m. 
Consequently, large uncertainties arise when picking 
traveltimes for greater separations. Furthermore, there 
is strong noise throughout the data for the third geophone 
of Apollo 14 ASE thumper shots. Strong reveberations exist 
in the data of thumper shots 1 through 7 for the second 
geophone of Apollo 14 ASE. The first breaks of the signal 
onsets can not be determined with any accuracy for these 
data. Also, the amplitudes of the first arrivals of Apollo 
16 ASE thumper shots are smaller than those of Apollo 14 's 
at the same geophone/source distances. 

There is little noise in the grenade launching data. 
However, the onsets on the third geophone (95 m) of all the 
grenade launchings are still too low to determine first- 
arrival times. The traveltimes of the first arrivals are 
expected to be the same at each peophone for these three 
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grenade launchings because they have the same geophone/ 
source distance to each geophone. But we find a maximum of 
5^ec difference in the traveltimes between Grenade 


— 

Launchings 2 and 4. The travelti.e differences might be 
caused by the difficulties in recording the initiations of 
the launchings. For the ASE grenade data, the first 


arrivals at all three geophones merge into the background 
noise so that the traveltimes of the first arrivals can 
not be determined. 


Stackin g and Anmi < . 

To further improve the data, those traces which have 
the same geophcne/source distances were stacked together. 
The stacking would be destructive to the noise and 
constructive to the signal because of the randomness of 
the noise from trace to trace and the coherence of the 
signal. Consequently, stacking the traces for equal 
separations will increase the amplitude of the signal by 
the number of traces stacked while the noise will increase 
by the square root of the same number. This gives a S/N 

ratio improvement equal to the square root of the number 
Of traces stacked. 

Per the ASE thumper-shot data, those traces which have 
the same geophone/source distances are listed in Table 10. 
The numbers listed on the right-hand side of the table are 
the thumper-shot numbers which have the given separations at 
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their left. However, those thumper shots with an aaterick 
above their numbers (ten are from the Apollo 14 ASE and one 
is from the Apollo 16’ s) have such high noise levels that 
they were not used in stacking. The number of traces 
stacked are given in the second column of Table 10. There 
are between 2 and 7 traces with the same separations for 
the ASE thumper shots. The stacking would give S/N 
improvements between ^2 and/?. The stacked data are then 
amplified so that they are plotted at almost full scale for 
each stacked trace. This helps us pick the first arrivals. 
The amplifications used are given in the third column of 
Table 10. 

Figure 22 shows the results of the stacked, filtered, 
and amplified profile of Apollo 14 and 16 ASE thumper shots. 
Remember that the thumper firing time is 1.21 seconds 
after the beginnings of the traces , and the increment of 
the geophone/source separation between traces is about 
4.57 m. The first arrivals could be determined easily 
and accurately for separations up to 32.0 m. 

Figure 22 also shows a smooth curve for the first 
arrivals on the profile which suggests that the continuous 
velocity variation model might be the better one. If the 
velocities increase in a stepwise fashion, the traveltime 
curve of the first arrivals on the profile will then be 
straight lines. However, it was still not possible to pick 
first-arrival times for distances greater than 32.0 m. 
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Figure 22. Stacked, filtered, and amplified 
profile of Apollo ASE’s, 
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The stacked profiles and the single geophone profiles 
of Apollo 14 and 16 AS£ thumper shots are compared with the 
results of the stacked profile of the Apollo 14 and 16 ASE 
thumper shots. If there are any significant differences 
in the velocity structures at these two sites, the results 
of these profiles will show the differences. Little 
difference is found. 

Traveltime Variations t 

Using Kaufman’s (1953) result, a 'velocity function 
with depth, z, given by V(z)sV^(z/z^)^ has a traveltime 
function, t(x), with separation, x, given by 
t(x) = cx"/Vq 

where c is a constant, m equals 1-n, and is the velocity 

at a depth of z^. The above equation also represents the 

travel time/separati on function of the direct waves for 

the homogeneous and layered model (m=l); where c/V =1/V 

0 

and V’ is the constant velocity of the surface layer. 

The slope of the traveltime data plotted on a log t 
(traveltime) versus log x (geophone/source separation) 
graph would determine the exponent, m. 

The travel time data of the first arrivals measured 
from the Apollo 14 and 16 ASE stacked profiles for 
separations up to 32.0 m are listed in Table 11 (from Gangi 
and Yen, 1979). The questionable data are given a weight 


Table 11 
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Of one quarter that of the hlgh-S/N data, the V -a are 
oaxculated from the slopes of the least-square-mted 
straight lines and extrapolated to a„=l 

V„ at a depth of i „oes not imply that the powdered 
layer would extend to 1 km. The at 1 km is. for 

convenience, a constant of reference used to characterize 
the velocity model. The depth, z„. can be chosen at any 
other depth, say z then the V^'s in Table 11 will be 
multiplied by Zj - . Table 11 (from Gang! and Yen, I979) 

Shows a greater variation of the reference velocities, 

Vo's, than of the slopes. The reference velocities 
vary between b30 m/sec and 630 m/sec while the slopes 
vary between 0.74 and 0.80. While the measured slopes are 
variable, they are all consistently lower than m«l-n«l, 
the value that would be obtained for the homogeneous and 
layered model. They tend to the value predicted by the 
self-compactlng-powder model, namely, m.l-l/6.0.833. 

The traveltimes of the first arrivals calculated 
for the Kovach and Watkins models (I973) at the Apollo 14 
and 16 landing sites are tabulated in Columns A and B of 
Table 11 (from Oangi and Yen. 1979). There are not large 
differences between the measured traveltimes and the 
calculated traveltimes. However, the biggest d.-fferences 
between the Kovach and Watkins models and the measured 
values occur at the smallest separations. They assumed 
that the traveltimes at zero separations are zeros, and 
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they fitted most of the data points simply by forcing the 
least-squared travel time curves to pass through the zero 
separations (Figure 2, page 5). 

In stacking, we assume that the separations are 
exactly the same for each trace. Stacking should 
improve the S/N ratio by the square root of the number of 
traces stacked. However, the S/N improvement obtained 
from the stacked data is not that expected theoretically. 
This might be due to the time offsets introduced into the 
individual traces by differences in separations and small 
differences in elevations at the source and receiver 
locations. Figure 23 shows the second geophone profile 
of the Apollo 14 ASE thumper data (the firing time is 
approximately 1.21 seconds after the beginnings of the 
traces). Thumper shots 14, I5, and 16 of Apollo 14 ASE 
(corresponding to separations of I3.7I m, 18.29 m, and 
22.86 m, respectively) misfired. Notice that some good 
first arrivals can be found at separations greater than 
32.0 m (Figure 23). However, strong reveberations exist 
in the data of thumper shots 1 through 7> the first 
arrivals could be detected but the accuracy of the onset 
determinations is very questionable. 

In Table 12, the traveltime data from the first and 
second geophone profiles of Apollo 14 ASE thumper shots 
up to 45.72 m are tabulated) that is, traveltimes from 
thumper shots 11 through 21 are listed except for those 
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Figure 23. Filtered and amplified profile of the 
second geophone of Apollo 14 ASE. 
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Table 12. 


thfn?Bt“arr“afs!“"®‘’ "''“''"“I”*- 


Of 


B 


Separation, m 


tt 


#* 


4.57 

9.14 

13.71 

18.29 

22.86 

27.43 

32.00 

36.58 

41.14 

45.72 


Slope 
, m/se( 


Traveltime, msec 


Measured' 


GP 1 


53 

94 

132 

164 


* 

# 

312 


GP 2 


53 

94 


230 

260 

276 

293 

315 


0.82''’ 

254‘'''’’ 


Calculated 


53 

95 

132 

168 

201 

234 

*259 

*277 

*295 

*313 


0.82 

254 


B 


52 

93 

131 

166 

200 

233 

*259 

♦277 

*295 

*313 


0.833 

254 


Measured from the first and second peophone 

profiles of Apollo 14 ASE . Data are filtered 
between 20 and 40 Hz. nj-tered 

For y-373(z/zo )^ * in the surface powder laver 
overlying a homoreneous medium with V.=254 m/sec. 

surface powder layer 

Misfire 
Noisy data 

Least-square fitted, for separations less than to 

3^0 Paratl ons between ^ 

j<:.o m and 45.72 m, inclusive. 

The thickness of the powder layer is 11 m. 


m 
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misfired. Beyond those separations, the first arrivals 
are ambiguous. The third-geophone profile of Apollo-14 
ASE thumper shots is not used because of the low S/N ratio. 
On a linear travel time/separation plot (Figure 24), the 
first arrivals show a smooth curve up to approximately 
30 m» no straight line could be fitted to those data 
points which would pass through the origin in the same 
time. It strongly suggests that the power-law velocity 
model might be a more suitable representative for the 
velocity structure on the very shallow lunar crust. 

Beyond 30 m, a straight line is observed. Those data 
points are assumed to be associated with refracted waves 
from the second layer which is assumed to be homogeneous 
with a constant velocity of 254 m/sec. However, the 
velocity in the second layer is still questionable because 
there only are five data points and uncertainties increase 
in the data at larger separations. 

The log/log traveltime plot for the single geophone 
profiles of Apollo 14 ASE thumper shots is shown in 
Figure 25 along with the least-square-fitted line with a 
slope of 0.82. It can be seen from the figure that the 
straight line is an excellent fit to the data. The 
reference velocity, V^=373 m/sec, is evaluated from the 
slope of the least-square-fitted line and the time 
intercept at x*l m. If the exponent of the travel time 
function, t(x)*cx®AQ (where x is in meters), is 0.82, 


TRAVELTIME (msec) 


400 


300 




/ 


100 - / 

/ 

/ 


m 


10 20 30 


SEPARATION (m) 


Theoretical and measured traveltiroe curves 
Measured data are obtained from the first 
and second geophone profiles of Apollo 14 
ASE thumper shots 11 through 21. Direct 
waves are in dashed line. Refracted waves 
are in dot-dashed line. Reflected waves 
are in solid line. 



TRAVELTIME (msec) 




62 


c equals 5.77. However, the slope of the least-square- 
fitted line is almost identical to the value, m= 1-1/6=0,833, 
predicted by the self-compacting-powder model that the 
one-sixth power velocity variation is very possible. 

Using Kaufman's (1953) result, a velocity function with 
depth, z, given by 
V(z) = 

has a traveltime function, t(x), with separation, x, 
given by 

t(x) = 1.2(15 Fz/8)^/^x5/Vv^. ( 2 ) 

A least-squares fit of this function to the measured 
traveltimes gives V^=345 m/sec. The theoretical traveltime 
curve of the direct waves with V^=345 m/sec fits the data 
points very well (Dashed linej Figure 24, page 60). 

If this is the velocity model for the very shallow 
lunar crust, we can derive the traveltime equations for 
the refracted and reflected waves. That is, we assume 
there is a one-sixth power velocity variation oi a self- 
compacting-powder layer on the lunar surface with a 
reference velocity of 345 m/sec at 1 km and a thickness 
of H overlying a homogeneous layer with a constant 
velocity of 254 m/sec. For the travel time equation, to, 
of the refracted waves, we have 

\ - x/Vi.{ViV4V^^)F(e^) (3) 

XI separation 


where 
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VqI reference velocity of the powder layer 
at km 

I constant velocity of the homogeneous layer 

P(®c^ = 36^-sin O^cos 0g(3+2sin^ 6 -8sin^ 6 ) 

c c 

e^j« critical angle; where sin 
Hj thickness of the powder layer 
t^« intercept time 

The thickness, H=ll ro, of the powder layer is calculated 
from the V^, , and t^ using the fixed-point iterative 

method (see, for example, Conte and de Boor, I972, page 44). 
The maximum separation for the direct waves is 
x= 5.89H=64.8 m (Gangi, 1972). 

The traveltime equation of the reflected waves must 
be expressed in a parametric form. That is, the travel time 
and separation relationships, t^(p) and x^(p), for the 
reflected waves are (in terms of the ray parameter, p, 
which is given by the Snell's law, p=sin 6(z)/V(z)) 

tj,(p) = (3Zo/2pV^^)(3a^sin"^ (b/a)-(3a^+2b^)bc) 
*r(p) = (Zo/4)d5a^sin"^ (b/a)-(15a^+10aV+8b^)bc) 

(4) 

where a s i/pv 

^ o 

b ' 

c = (a^-b^)i 

The theoretical traveltimee for the refracted and reflected 
waves, dot-dashed and solid lines, respectively, are also 
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/[ self -compacting powder layer 

y\ V-345(z/l0OO)''* m/sec. ^ 


DEPTH* llm\ 


z is in m. 


DEPTH: UNKNOWN 


homogeneous layer 
Vi *254 m/sec. 


Figure 26. Velocity structure at the Apollo 14 landing 
site. 


shovm in Figure 24 (page 60) based on the values found. 
Figure 26 shows the velocity model at the Apollo 14 
landing site. Notice that the velocity at the lunar 
surface, z=0 m, is zero based on the velocity function, 
and the thickness of the second layer, the homogeneous 
layer, is unknown. However, whether the second layer 
is homogeneous is still unknown. 

For the first geophone profile of Apollo 16 ASE, the 
first arrivals can be picked for separations up to 18.29 m 
and, with difficulty, for 2?. 43 m (see Table 13). 

After 27.43 m, the first arrivals are buried in the noise. 
The slopes (see Table I3) of the least-square-fitted lines 
are, again, lower than that predicted by the Kovach and 
Watkins model (1973) • B'or the three single -geophone 
profiles, the variation of the reference velocities is 


raoie ij. 


iraveltimes measured from single geophone 
profiles of Apollo 16 ASE thumper shots. 
Data are filtered between 20 and 40 Hz. 
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Thumper shots 12 through I 9 
Thumper shots 1 through 11 
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greater than that of the slopes i the values of the 
reference velocities vary from 1059 m/sec to 773 m/sec. 

The values of the reference velocities are almost doubled 
or tripled that of Apollo 14 ASE. The slopes of the 
least-square -fitted line vary from O.65 to O.71. 

Empirically speaking, the velocity variation on the lunar 
surface at the Apollo 16 landing site tends to the fourth- 
root velocity variation rather than the sixth-root velocity 
variation. However, if we also assume that the velocity 
variation on the lunar surface at the Apollo 16 landing 
site is one-sixth power, a least-squares fit of the 
travel time function (Equation (2), page 62) for the direct 
waves to the measured traveltimes gives m/sec. 

The traveltimes for the Apollo 16 ASE grenade 
launchings are tabulated in Table 14. Notice that there 
is a maximum time offset of 57 msec between Grenade 
Launchings 2 and 4. The time offsets between the other 
two pairs of grenade launchings are approximately 30 msec. 
The same traveltimes are expected at each geophone 
because the geophone/source separations are the same. 
Furthermore, the first arrivals for Grenade Launchings 3 
and 4 at Geophone 3 are lu-.de tec table , and the first 
arrival at Geophone 3 of Grenade Launching 2 is also 
questionable because of the extremely low amplitudes. 

Based on the previous knowledge obtained, we assume that 
the direct waves are the first arrivals at Geophone 1 (14 m) 
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Table 14. 


Travel tiroes measured from grenade 
launchings of Apollo 16 ASE. Data 
are filtered between 3 and 4o Hz. 



while the refracted waves are the first arrivals at 
Geophones 2 (50 m) and 3 (95 m). We find the velocity, 

302 m/sec, of the refracted waves from the travel time data 

of the second and third geophones of the second grenade 
launching. 

An attempt was made to correlate the grenade launching 
data with the results of the thumper-shot data of 
Apollo 16 ASE. Using the reference velocity, 357 m/sec, of 
the powder layer and holding the traveltime differences 
among those at the geophones of the grenade launching to be 
the same, we correct the travel time at Geophone 1 for 
Grenade Launching 2 to 129 msec, determine the intercept 
time of the refracted waves to be 145 msec and find the 
thickness, H*12 m, of the (surface) powder layer. 

It is of some interest to notice that the deviations of the 
velocity structures between the Apollo 14 and 16 landing 
sites are within 16^i namely, 4?5 for the reference 



velocity, 8^ for the thickness of the powder layer, and 
16^ for the constant velocity of the homogeneous layer. 
However, it is not surprising that a larger deviation occur 
for the constant velocity of the homogeneous layer because 
we essentially have only two data points to be interpreted 
for the constant velocity of the refracted waves, and the 
one at Geophone 3 of Apollo 16 ASE grenade launching is 
very questionable because of the extremely low amplitude 
of the first arrivals. If we let the constant velocity, 
of the homogeneous layer be either 302 m/sec or 254 m/sec 
and vary the reference velocity, V^, of the powder layer 
from 300 to 420 m/sec, we find the depth varies from 11 
to 13 m and 9 to 10 m, respectively (Table I5). 

Amplitude Variations » 


The amplitudes of the first arrivals for separations 
up to 32.0 m are measured from the profiles. The 
measurements of the amplitudes for the ASE data are more 
difficult than that of the travel times. These difficulties 
(summarized by Gangi and Yen, 1979) are caused byi 1) the 
coarseness of the amplitude sampling, 2) the variabilities 
of the thumper-shot strengths, 3) the variabilities of 
the geophone siting and coupling, and 4) the low S/N ratio 
for larger separations. However, the interpolating effect 
of the Butterworth bandpass filter reduces the coarseness 
of the amplitude data, the stacking reduces the 
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Table 15. Calculated values of the Intercept time, 
Tq. and their related depth, H. 


*1 

'^C- 

... / 

V^*302 

tn/sec 

Vj^» 254 in/sec 

msec 

m/sec 

T^, msec 

H, m 

Tq, msec 

H, m 

150 

139 

128 

118 

110 

300 
1 330 

360 
390 
420 

169 

155 

144 

134 

126 

11 

11 

12 

12 

13 

137 

123 

112 

102 

94 

9 

9 

9 

10 

10 


^ ’ i tHight® 1 l>ased on the on 



variabilities of the shot strengths sna the geophone 

aensitivitles. The amplitude data, sgain. are plotted on 

a log/log graph and compared with the amplltude/separatlor 
variation predicted by the two models. 

Since the thumper shots give primarily vertical 
forces and the geophones are vertically oriented, the 
amplitude function of the direct arrivals for the 
homogeneous and layered model is given by 
Mx) . A^(x/x^)-2 

for small separation (see. for example. White. I965. 
page 215). while the amplitude/separation variation for 

the power-law-velocity model is estimated to be (Gang! and 
Yen, 1979) 

A(x) . A„(x/x„)-(13-a)/12 , 3^1 

where s is a measure of the source radiation pattern in 

the vertically inhomogeneous medium. The decrease of the 

amplitude with separation of the direct waves in the 

powder-layer model is less than that in the homogeneous 
and layered model. 

The measured amplitudes along with the slopes of the 
least-s,uare-fitted lines are given in Table 16 (from Gangi 
and Yen. 1979). goth the amplitudes for the single geophone 
profiles and the stacked profiles are given in the table. 
Measurements are made on the data that have been bandpass 
filtered with -3 dB frequencies of 3 to 66 Ht and 20 to 
40 Hz. Notice that there is a great deal of scatter in the 


Table 16. 

A. aandputcd; S-MHt 


(arbitrary units 
(From Ganpi and Yen, I979) 
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(^ta. Some of it ie becauee of the variabilities of the 
shot* strengths and geophone sitings, while some is due 
to the low S/n ratio and the coarseness of the amplitude 
sampling. Figure 2? (from Gangi and Yen, I979) ahows a 
representative amplitude plot (for Geophone 3, Apollo 16 
ASE). While all the data are well fitted by the straight 
line with a slope of -2.01 in the figure, the two largest 
amplitude data, which have the highest S/N ratio, suggest 
a lower slope. 

Table 16 (from Gangi and Yen, 1979) shows the slopes 
for the fairly good data lie between -I.5 and -2.5. 

The slopes are closer to the slope predicted by the 
homogeneous and layered model than that predicted by the 
powder-layer model. However, the slopes predicted by 
these two models are made on the following assumptions t 
1) all the thumper sources are of equal strength, 2) all 
the geophones are equally coupled and oriented to the 
lunar surface, 3) there is no attenuation, 4) there is no 
energy loss by conversion of P-wave energy into S-wave 
energy, and 5) there are no scatterers in the lunar 
regolith. The first two effects would increase the 
scatter in the data, while the latter three effects would 
decrease the amplitudes. Therefore, the two predicted 
slopes should be considered as upper bounds on the measured 
data. However, energy loss does exist when seismic waves 
propagate through a medium. The energy loss is called 
the attenuation or absorption. In general, the attenuation 


AMPLITUDE (ARBITRARY UNITS) 


6E0PH0NE 3 
APOLLO- 16 ASE 
BANDPASS* 3-66 Hz 
SLOPE* -2.01 


I 5 6 7 8 9 10 20 3C 

SOURCE /GEOPHONE SEPARATION (m.) 

Figure 27. Log/log plot of the amplitudes 
versus separations. 

(From Ganpi and Yen, 1979) 
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in a medium is an exponential function with distance 
(see, for example, Dobrin, 1976, page 39) as follows 
A(x) = A^x"*^e"^* 

where Q is the coefficient of attenuation. Using the 
amplitude data cf the third geophone profile of Apollo 16 
ASE thumper shots, we find (also see Figure 28) 

A^ = 71.66 
n = 1.463 
a = 0.04? m“^ 

For the powder-layer model, the exponential factor, -n, 
results in a negative value for s which contradicts the 
condition s>l (also see Equation ( 5 ), page 70 ). However, 
the exponential factor, -n, is close to but less negative 
than the slopes shown in Table 16, and the predicted 
values for the powder-layer model and the homogeneous and 
layered model are treated as upper bounds on the measured 
data. This value, -n=-1.463, seems to favor the powder- 
layer model. Figure 28 shows that the amplitude decreases 
exponentially at larger separations, while the factor 
x dominates the amplitude variation at small separations. 

Velocity-Spectrum Analyses i 

The velocity-spectrum technique used in this research 
was pioneered by Taner and Koehler (I 969 ). The semblance 
of the velocity spectrum shows the power of the signals 
arriving at different vertical two-way traveltime, t^, over 
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A (X) ■ AoX*"e*®*^ 



Figure <.8 . Amplitude attenuation with separation. 
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a spectrum of (assumed) root-mean-squared (RMS) velocities, 
V. This technique is used to search for the parameters of 
the (assumed) velocity model for the direct, reflected, 
and refracted waves by delaying and summing the traces 
along a profile. The delay is computed on the basis of 
the (assumed) velocity model and the maximum amplitude 
for the sum over a time window is obtained when a proper 
velocity model is used. That is. the semblance is 
calculated for the digitized data by 


S(t^.V) = 


M 


Jq+N M 

(Za..)' 


Jo^N 

2 


M 


Za, ,2 

i=l 


( 6 ) 


where Is the amplitude of the J-th sample point of the 
i-th trace on the profile, M is the number of the traces, 
and N is the width of the time window. The sampling 
point, J^, in the i-th trace is obtained by 

= INT( T(t^,V,X^)/At ) (7) 

where t^^ is the vertical two-way traveltime or the 
intercept time, V is the RMS velocity, is the geophone/ 
source separation for the i-th trace. At*1.89 msec/sample 
is the sampling interval, INT represents INTEGER in 
computer language, and T is the traveltime function based 
on the (assumed) velocity model and the event of interest. 

The velocity spectrum is generated by measuring the 
semblance, using Equation (6), for various t^ and V. 


calculated for v 

v...tov . Of .7) 
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to V^,. intercept tine is thenTnc^asril 

k 4> ^ ^ e 


4 increasArf 

ncremente of 4t„ and the same procedure is repeated 

The maximum semblance value will occur at the v , ’ 

to and V associated with the traveUlme ^ 

event (Figure 20 . the 
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or 

Notice that the maximum value of the semblance is 1 when 
there is only signal in the time window and its minimum 
value is 0 when there is no signal but noise. The S/N 
ratio can then be estimated by taking the square root of 

V*’n- ’'•'O semblance read on the peaks of the velocity 

spectra varies from approximately 0.2 to 0.7 which 
correspond to S/N ratio of 0 . 5 to 1 . 5, respectively. 

The assumed velocity model based on the previous 
investigation is that there is a self -compacting, powder 
layer overlying a homogeneous medium. For the direct 
waves, the traveltime function T(t^,v,xp of Equation (7) 

(page 76) has the form 

T(t^,V,X^) = AX.^/V V 

where t^ equals zero. A is a constant (also see 
Equation (2). page 62). and V. in this case, is the 
reference velocity at 1 km for the powder layer. Also, 

Gangi (I972) demonstrated that the travel time, T^, for 

a ray reflected n times (the n‘‘=^^-surface-reflected waves) 
from the surface is 

Where T^ is the travel time of the primary direct waves. 

The above equation shows that the apparent velocity of the n- 
th-surfaoe-refleoted waves is equal to the product of 
( 1 +n) and V for the direct waves. 

For the reflected waves. Taner and Koehler (I969) 
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..v> 








IT IT -^-tion (.,. 

P ge 63) can be approximated by 


where 


*t(x) « (c^+C2X^+c.x^+ 
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®1 * 

°2 “ ®l/®2 “ ^ 


^ = v-2 

C 3 = (a2^-aj^a^)/4a2^ 
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i 


= 2 
m 


f" v 2 "- 3 (,m, 

V(z), the sixth-power velocity function 
with depth 

Hi the depth of the powder layer 
I^«•the^nore. they conclude that the first-two-ter™ 
apprcxination of the above e,uation will converge ,uite 

ITeT ^'•aveltin.e/separation 

curve for most cases of T>ra+if»oi • 4. 

5 01 pratical interest. The first- 

two- term approximation has tha «=i«xa 4. 

nas tfte same form as Dix’s 

formula. That is, 

T(t„,V,X^) . (t^ 2 .x^ 2 / ^ 2 )i 

For (n-i)-th multiply reflected waves the eh 

waves, the above equation 

is also valid. Its intercept time, t is „ « 

the primary reflected waves, while the“’ ‘ 

remains the same. The relationship 

ne relationships among the c. *s, t 

and V can be obtained immediately as follows ^ 


c i - + ^ fH 

1 ■ ^o.n ° 2nJ dzA(t) 
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or 


«2 ‘ V(i)dz/ 2nj"" <Jz/v(8) 

2n fH 

“TTJ V(z)dz 

0,1 ■* 0 

= n( 72/245 ^ 


:4.210V^/^(t 

0 I 1' 


H ^O.lf93Vt„_j/n 
'^0 

"here to.n*"\.l- 0'>''i<»»sly, we would expect to have the 

semblanoe peaks at every grid point of (nt^_j.7') for 

the (n-l)-th multiply reflected waves, where t and V- 
are the Intercept time and the (assumed) RMS velocity. 

respectively, of the primary refiected waves and n is 
equal to or greater than 2. 

For the refracted waves, the traveltime function 
Ktg.V.Xj^) (also see Equation (3), page 62) is 
T(t„.V.X^) = t„^xj/ V 

Where V is the constant velocity, Vj. of the homogeneous 

medium. This verocity remains unchanged for the n-th 

multiply refracted waves, while its intercept time is 

n-n times that of the primary refracted waves. However. 

the intercept time is a transcendental function in terms 

of V^, H, and Vj. Either or H has to be assumed to 
determine the other one. 

A synthetic seismogram and computer program have 
been used to test this velocity-spectrum technique. 


The synthetic seismogram consists of the primary direct, 

surface reflected, reflected, and refracted waves t the 
assumed velocity model is a powder layer, with a thickness 
of 10 m and a reference velocity of 330 m/sec at a depth 
of 1 km, overlying a homogeneous layer with a constant 
velocity of 250 m/sec. The geophone/source distances of 
the synthetic profile are from 4.57 m, in increments of 
4.57 m, to 45.7 m. The arrival signals are one and half 
cycles of a sine function with a duration of approximately 
57 msec (Figure 30 (a)). Random noise is introduced in 
the synthetic profile. Figure 30(b) shows a synthetic 
velocity spectrum of the direct waves. The first-sur- 
face-reflected waves are contained in the profile and 
its reflection coefficient is 0.35* The S/N ratio is 
approximately 2.6. The solid, dot, and dash curves 
represent the semblances for time windows with durations 
of approximately 19, 38, and 57 msec, respectively. The 
spectrum indicates that the longer the duration of the 
time window, the lower the semblance value at the peak 
and the wider the spread of the semblance peak. Notice 
that the semblance curves shift to the right at the 
neighborhood of the velocity, 330 m/sec, of the primary 
direct waves and shift to the left at the velocity, 

294 m/eect of the surface-reflected waves as the duration 
of the time window increases. These tests also show 
that interfering signals on one or more seismic traces or 
high noise may cause spurious peaks and lower the 
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Figure 30. 


VELOCITY, m/sec 
(b) 

arrival signal which is one and 
half cycles of a sine function with a 
approximately 57 msec, 
kb; The velocity spectrum of the direct 
waves for the synthetic profile. The 
solid, dot and dash curves represent the 
semblances for time windows with duratio 
of approximately I9, 3^, and 57 msec, 
respectively. 
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semblance values at the peaks on the velocity spectra 

This would explain the low semblance values and the 

spurious peaks on the velocity spectra of the Apollo 
ASE data* 


Figure 31 shows a representative velocity spectrum of 
the direct waves for the ASE data which is from the 
second geophones profile of Apollo 14 ASE. Since the 
semblance is very sensitive to the interference and 
noise, only clear traces can be usedi that is, thumper 
shots 12 through 21 are delayed and stacked (also see 
Figure 23, page 57). The solid, dot. and dash curves, 
again, represent the semblance curves for the time 
windows with durations of approximately I9, 38, and 
57 msec, respectively. The semblance peaks near 
370 m/sec are believe to be spurious because the peak 
of 19 msec time window is not stationary for the other 
two time windows and the semblances drop almost by a 
factor of 2 in that region. However, the semblance 
curves at the neighborhoods of approximately 320 and 
282 m/sec show similar character to that found in the 
synthetic velocity spectrumj that is, the curves 


shift to the right and left at the neighborhoods of the 
velocities of the primary direct waves and the surface- 
reflected waves, respectively, when the duration of the 

time window increases (Figure 30 (b)). This similar character 

suggests that the values of 320 and pfto «/« 

and 282 m/sec may correspond 

to the velocities of the primary direct and the first sur- 




Hce reflected waves • respectively. Furthermore, if the 
value of 320 m/sec is the velocity of the primary direct 
waves, the velocity deviation between the theoretically 
expected value for the first surface-reflected waves and 
282 m/sec is only about 1 %. However, the determination 
01 the velocities for higher orders other than that of the 
first surface -reflected waves is difficult because their 
apparent velocities are so close to each other that the 
determination from the velocity spectrum is impossible. 

Figure 32(a) and 32(b) show the synthetic velocity 
spectra, using a I9 msec time window, for the reflection 
and refraction, respectively. The expected values of 
the intercept times and the RMS velocities for the 
synthetic reflected and refracted waves are I57 msec 
and 130 m/sec, I34 msec and 250 m/sec, respectively. 

There are offsets between the locations of the expected 
peaks, marked "o", and the measured peaks, marked 
found in the synthetic velocity spectra. The offsets of 
peak locations are caused either by the interfering of 
the seismic signals or the noise. However, the centers 
of the semblance contour are not far away from the 
expected ones; therefore, the semblance -contour centers 
for the reflected and refracted waves of Apollo ASE data 

will be used to evaluate the parameters of the assumed 
velocity model. 

Figure 33 shows a representative velocity spectra 
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Synthetic velocity spectrum for the 
reflected waves using a 19 msec time 
window. The S/N ratio is 2.6. The 
dot represents the measured semblanc 
peak and o represents the expected 
semblance peak. 
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Figure 33. The velocity spectrum for the reflected 

geophone profile 

of Apollo 14 ASE using a I 9 msec time 
window. The dots indicate the semblance 
peaks and x represents the center of the 
semblance contour. 
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for the reflected waves of the Apollo ASE data. The 
velocity spectrum is from the first geophone profile of 
Apollo 14 ASE using a time window of 19 msec duration. 

The dots indicate the semblance peaks and x indicates the 
center of the semblance contour. The contour shows a 
similar charaoter to that of the synthetic profile. The 
time windows of longer duration will only broaden the 
contour of the semblance peak but have similar character. 
The intercept time and the RMS velocity at the contour 
center of the velocity spectrum are at approximately 
140 msec and 121 m/sec. respectively, which gives for 
the thickness. H. and the velocity. V„. 8.4 m and 
318 m/sec. respectively. This independent evaluation 
verifies-4;he reference velocity of the surface powder 
layer. V^. is approximately 320 m/seo. 

Figure y, shows the velocity spectrum for the 
refracted waves which is from the second geophone profile 
Of Apollo 14 ASE using a time window of 19 msec duration. 
The dots and x. again, represent the semblance peaks 
and the center of the semblance contour, respectively. 

There are only seven clear seismic traces on the second 
geophone profile of Apollo 14 ASE (also see Figure 23. 
page 57). However, those traces with geophone/source 
distances less than the critical distance can not be 
used for the velocity spectrum analysis of the refracted 
waves because refracted waves do not exist before the 
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critical distance. The critical distance is assumed to 
be less than I3.7I m. That is, there are only five 
available traces for this analysis of the refracted waves 
The semblance peak is distributed in a very narrow belt 
area. The region for the intercept time and refraction 
velocity is from 120 msec and 213 m/sec to I35 msec and 
246 m/sec, respectively. If we assume the reference 
velocity of the powder layer be 320 m/sec, then the 
corresponding thickness of the surface powder layer vary 
from 9.0 to 9*7 m, respectively. Furthermore, the 
center of the semblance contour for the refracted waves 
is approximately at 128 msec and 230 m/sec which 

corresponds to a thickness of the surface powder layer 
be 9.3 m. 

The previous investigation in the section of 
"Traveltime Variations” show that the-reference velocity 
of the surface powder layer is approximately 350 m/sec. 

The thickness of the powder layer is approximately 11 m. 
The constant velocity of the homogeneous medium is 
250 m/sec. The results of the velocity-spectrum analysis 
are consistent with those values found in the previous 
investigations. The differences of the reference 
velocitjr, the thickness of the powder layer, and the 
constant velocity of the homogeneous medium are 9 » 6 %, 

15 *5^* and 8.0?5, respectively. These percentages of the 
differences are reasonable and possible. 



G e^hone-Cq a pj.infi ^n4_Shot-Strength Variah< 


93 


An analysis of the amplitTide data is made to 
investigate the exponent of the amplitude variation and 
the variabilities of the geophone sensitivities (in place) 
and the thumper-shot strengths, if there are sufficient 
redundancy in the data. The amplitude variation of the 
direct waves can be written in a general form both for 
the constant velocity medium and for a powder-layer 
medium as follows (Gangi and Yen, 1979) 


'ij “ Vj 


where 


G^t the sensitivity of the i-th geophone at 


Sjj the strength of the j-th thumper shot at X. 
m » the exponent of the amplitude variation 
The above equation is normalized and linearized in terms 
of the relative geophone sensitivities (g^=G^/Gj), 
relative thumper-shot strength (s.=S./Sj), the exponent, m, 
of the amplitude variation, and an arbitrary constant by 
taking the logarithm. The arbitrary constant is the product 
of the reference geophone sencitivxty and the reference 


thumper-shot strength (a^«GjSj). Minimizing the summed, 
weighted, and squared error function as a function of 
those parameters (g^'s, s^’s, a^, and m), the above 
equation results in a matrix form (for details* see 


Gangi and Yen, 1979) 



( 8 ) 


9 ^ 


• • • » 


where a is a vector, whose components depend on the measured 
Mplltudes. the weights, and the geophone/souroe distances, 

A is a squared and symmetric matrix whose components 
depend on the weights and distances while p is a vector 
Whose components are the unknown parameters, gj .... a 

®o* assume that the distances are measured ■ 

with high accuracy and we establish the weights which 
depend on the qualities of the date. Consequently, the 
stability of the inverse matrix and the errors of 
the vector 5 are established only on the accuracy of 
the amplitude measurements. 

The weights for the Apollo 16 ASE thumper shots 
are tabulated in Table 17 (from Gangi and Yen. I979). 

Table 17 (from Gangi and Yen, 1979) shows that 19 of I9 
thumper shots give useful amplitude data and only 3 
(Thumper-shots 6. 7. and I7) of these 19 thumper-shots 
give first-arrival amplitude data on more than one geophone 
(also see Table 16, page 71, from Gang! and Yen, I979). 
Geophone a and Thumper-shot 7 are used as the reference 
geophone and the reference thumper-shot, respectively. 
Consequently, we have six amplitude measurements (two for 
each thumper-shot) to solve six unknowns a,, gj, g . s., 

Sjy. and m. The matrix equation (Equation (8)) - in terms 

of the weights, distances, and correlated and measured 
amplitudes - becomes 


Shot No . , j 


Geophone No., i 








96 


“-3.734“ 


"3.50 

1.00 

.75 

1.00 

1.25 

10.35“ 

-0.211 


- 

1.00 

0 

0 

1.00 

2.62 

-1.164 


- 


.75 

.50 

0 

2.39 

-1.592 


- 

- 

- 

1.00 

0 

3.13 

-0.914 


- 

- 

- 

- 

1.25 

3.48 

2ll.710_ 







30.84_ 









Solving this rfiatrix equation, we find that the 
relative geophone sensitivities and thumper-shot 


strengths are 



gl 

^3 

®17 


Gj/Gg = .724 

G^/G2 ® 1.40 

s^/s^ = .803 
= .848 
m = -3.57 

Unfortunately, those values seem to be unreliable. The 
reason for the unreliability of the values is because of 
the lack of redundancy and the quality of the data. 

The 30 to hOjS differences in the relative geophone 
sensitivities are higher than that expected, but are not 
too unreasonable. The 15 to 20J5 variations of the relative 
thumper-shot strengths are possible but larger than 
expected. The values of the exponent, m. is almost 
doubled compared to that of the single geophone profiles 

and stacked profiles (also see Table 16 . page 71 1 from 
Gangi and Yen, I979) . 

A similar treatment can not be made on the 
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amplitude data of Apollo 14 ASE thumper shots because 
there are no correlated amplitude data for Geophones 2 
and 3 of Apollo 14 ASE due to the misfires and poor 
S/n ratio. Neither can it be performed on the grenade 
or grenade launching data of Apollo 16 ASE due to 
the difficulties in measuring their amplitudes. 
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CHAPTER V. SUMMARY AND CONCLUSIONS 


The seismic data used to determine the velocity 
structure of the very shallow lunar crust are from 
thirty-two thumper shots (thirteen from the Apollo 14 
ASE and nineteen from Apollo 16 's) and three grenade 
launchings of Apollo 16 ASE. The data are used to 
compare two velocity models; namely, the homogeneous 
and layered model (Kovach et al , 1971, I 972 , 1974) and 
the self-compacting-powder layer model (Gangi, I 972 ). 

To cover the maximum dynamic range, the seismic 
data are log compressed into thirty-two binary levels 
for the ASE's which gives a coarseness to the amplitude 
sampling. Furthermore, severe glitches are found in the 
lunar seismograms (Figure 16, page 39 ). They are 
recognized by the fact that they are of short duration 
and have values which are inconsistent with the preceding 
and/or following sampling values. 

To improve the quality of the data, they are 
"deglitched" to remove the extraneous values and filtered 
by four-pole, anti -aliased Butterworth bandpass filters 
with -3 dB frequencies at 3 and 66 Hz, and 20 and 4o Hz 
(Figure 21 , page 47). In addition, the Apollo 14 and 16 
stacked profile, the Apollo 14 stacked and Apollo 16 
stacked profiles, and the single geophone profiles of 
Apollo 14 and 16 ASE's are used to examine the velocity 
variation on the lunar crust at these two landing sites. 
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If there are any significant difference in the velocity 
structures at these two sites, the results of these 
profiles will show the differences. Little difference 
is found. 

For the traveltime variation, Gangi and Yen (1979) 
indicate that the travel time data plotted on the log/log 
graph will give the exponent of the velocity function. 

That is, the travel time/separati on function, in general, 
has the form 

t(x) = cx"*/V^ 
o 

where c is a constant, is the reference velocity at 
1 km, and m=0.833 predicted by the self-compacting 
powder layer or c/V^=V*, where V is the constant velocity 
on the surface layer, and m*l predicted by the homogeneous 
and layered medium. 

The traveltimes of the first arrivals (Table 11, 
page 5^5 modified from Gangi and Yen, 1979) can be 
determined accurately only up to separations of 32.0 m 
for the stacked profiles. While the reference velocities 
vary between 430 nj/sec and 630 m/sec for the stacked 
profiles, the values of the exponent, m, vary between 
0.74 and 0.80. The values of the exponent are all 
consistently lower than that predicted by the 
homogeneous and layered model and tend to that predicted 
by the powder-layered model. 

However, there are some good first arrivals 
found in the single geophone profiles at larger separations. 








. V. ■' 
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The traveltime data of the single geophone profiles of 
Apollo 14 ASE (Figure 24. page 60) show a smooth curve 
up to approximately 30 mt no straight line can be fitted 
to these data points which will pass through the origin 
at the same time. Furthermore, the value of the 
exponent, m=0.82, for those data points is close to that 
proposed by Gangi (1972). These strongly suggest that the 
continuous velocity model might be a more proper 
representative for the velocity variation on the lunar 
surface. We find the reference velocity at a depth of 
1 km for the sixth-root velocity variation to be 345 m/sec 
with a thickness of 11 m. Beyond 32 m, a straight line 
is observed and those data points are assumed to be 
refracted waves from the (second) homogeneous layer with 
a constant velocity of 254 m/sec (Figure 26, page 64). 

Unfortunately, the first arrivals for the Apollo 16 
single geophone profiles can not be determined with any 
accuracy for separations greater than 30 m because of 
the quality of the data. The values of the reference 
velocities vary from I059 m/sec to 773 m/sec for 
separations less than 30 ro, while the values of the 
exponent, m, vary from O.65 to 0.7I (Table I3, page 65). 

The exponent, again, is closer to that predicted by the 
self-compacting powder layer (Gangi, I972) than that 
predicted by the homogeneous and layered model (Kovach et al, 
1971, 1972, 1974). If we assume that the velocity of the 
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surface layer at the Apollo 16 landing site has a one- 
sixth power variation and we correlate the grenade 
launching data with the thumper-shot data of Apollo 16 
ASE, we find a reference velocity of 357 m/sec at 1 km 
with a thickness of 12 m for the (surface) powder layer 
and a constant velocity of 302 m/sec for the (second) 
homogeneous layers It is of some interest to notice that 
the deviations of the velocity structures between the 
Apollo 14 and 16 landing sites are within 16^j namely, 

4^ for the reference velocity, 8?S for the thickness of 
the (surface) powder layer, and 16fi for the constant 
velocity of the (second) homogeneous layer. 

For the amplitude variation, Gangi and Yen (1979) 
indicate that the theoretical amplitude/separation 
function for the homogeneous and layered model is x"^ 
while the approximated amplitude/separation variation for 
the self-compacting powder layer is where 

s>l, if we assume that* 1) there is no energy loss 
either by conversion or by attenuation, 2) no scatterers 
in the lunar regolith, 3) all the thumper shots are of 
equal strengths, and 4) all the geophones are equally 
coupled and oriented. The amplitudes of the firs s 
arrivals (Table 16, page 71 i from Gangi and Yen, I979) 
are measured up to separations of 32.0 m for all the 
stacked and single geophone profiles. However, the 
measurements of the first-arrival amplitudes are more 
difficult than those of the traveltimes because of the 


coarseness of the amplitude sampling and the poor S/N 
ratio at larger separations. The values of the exponent 
from the measured data (Table 16, page ?1, from Gangi 
and Yen, I979) all lie between -I.5 and -2.5. They 
tend to the value predicted by the homogeneous and layered 
model, not that predicted by the powder layered model. 
However ^he two predicted values of the exponent should 
be treated as the upper limits for the measured data. 

We-also find the value of the coefficient of 
attenuation,. 0.04? m“^, which dominates the amplitude 
variation at larger separations. The investigation of the 
variabilities of the geophone sensitivities and thumper- 
shot strengths (Gangi and Yen, I979) indicates that there 
are 30 to differences in the relative geophone 
sensitivities and I5 to 20^ variations of the thumper- 
shot strengths for the Apollo 16 ASE. These variabilities 
seem high but are not too unreasonable. This analysis 
requires accurate amplitude measurements and correlative 
amplitude data at two or more geophones. The same 
analysis can not be performed on the amplitude data of 

Apollo 14 ASE thumper shots because of the misfires and 
the poor S/N ratio. 

The semblances are calculated by delaying and summing 
the traces along a profile over a spectrum of the 
velocity, V, and the intercept time, t^. The maximum 
value of the semblance will occur at the values of t^^ and 
V Which are associated with the traveltime curve of the 
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event for the assumed velocity structure. The semblance 
is also related to the S/N ratio in the data. The maximum 
semblance is 1 when there are only signals. The semblance 
is 0, if the S/N ratio is zero. 

The semblances of the Apollo ASE data vary 
between approximately ©.2 and 0.? which suggests that 
the-S/N ratio is between approximately 0.5 and 1.5. 

However, spurious peaks are distributed over the velocity 
spectra. Tests of the velocity spectra of:. the synthetic 
seismogram indicate that the narrower the duration of 

the window t the higher the semblance value . Eurthermore 

the semblances of the velocity spectra are sensitive to 
interference of the signals and to noise; also, the 
semblance peaks shift to the right and left at the 
neighborhoods of the velocities of the direct and multiply 
direct waves, respectively. There ate offsets between 
the locations of the expected and measured peaks of the 
velocity spectra of the reflected and refracted waves. 

The centers of the semblance contours for the reflected 
and refracted waves are used to evaluate the parameters 
of the assumed velocity model. 

The assumed velocity model is that there is a powder 
layer whose velocity function is V(z)=V^(z/z^)^/^, where 
Vq is a reference velocity at a depth of z^, overlying 
a homogeneous layer with a constant velocity. The 
reference velocity at a depth of 1 km is approximately 


320 m/sec which is verified byt 1) comparing the 
velocity spectra of the direct waves with—that of the 
synthetic velocity spectra, 2) the velocity, 282 m/sec, 
of the first multiply direct waves which should be 2“^/^ 
times the reference velocity, and 3) the independent 
determim.tion of the reference velocity, 3I8 m/sec, 
from the velocity spectra of the reflected waves. 

The thickness of the surface powder layer evaluated from 
the velocity spectra of the reflected and refracted waves 
are 8.4 and 9-3 m, respectively. The constant velocity 
of the lower layer is 230 m/sec. These parameters of the 
velocity model are consistent with, those found in the 
previous investigation. All of them are within 10 %. 

An uneven powder-layer surface may explain some of the 
unusual characteristics found on the Passive Seismic 
Experiment seismograms (summarised by Gangi, I972). 

The one-sixth power velocity variation of the powder 
layer predicts that the velocity at the near lunar surface 
IS zero and the seismic rays return back to the surface 
at (or near to) 0° from the vertical. This will explain 
the lack of correlation between the vertical and 
horizontal components of the Passive Seismic Experiment 
seismograms. The uneven lunar surface of the powder 
layer will cause the "random walking" ol the seismic 
rays which, in turn, has a long duration of the signal 
(Gangi, I972i Gold and Soter, I970), 
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